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ABSTRACT 


^Plasticity  induced  closure  and  its  effect  on  the  fatigue  growth  of 
short  cracks  were  investigated  analytically  in  a  high  strength  titanium 
alloy  at  room  temperature  as  well  as  in  a  nickel  base  alloy  at  elevated 
temperature.  The  analysis  consisted  of  subjecting  a  single-edge  cracked 
specimen  with  an  initial  crack  length  of  .001  inch  to  cyclic  loading  and 
allowing  the  crack  to  propagate.  This  was  accomplished  using  a 
two-dimensional  finite  element  code  and  a  theoretical  finite  element 
model  of  a  typical  test  geometry. 

Visco  II,  the  two-dimensional  plane  stress/plane  strain  finite 
element  code,  uses  constant  strain  triangular  elements  and  incorporates 
the  Bodner-Partom  Viscoplastic  Flow  Law  to  handle  non-linear  material 
behavior.  This  is  a  time -dependent  material  model  which  can  easily  be 
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modified  to  handle  time- independent  behavior.  The  relations  in  this 
model  are  integrated  through  time  by  an  explicit  Euler  extrapolation 
scheme.  The  overall  solution  technique  involves  the  residual  force 
method.  Changing  boundary  conditions  are  incorporated  through  the  use 
of  a  crack  closure  algorithm  in  addition  to  a  crack  growth  procedure. 

'■—-'The  numerical  simulations  involved  subjecting  specimens  of  Tl-6246 
at  room  temperature  ( t ime - independent  behavior)  to  cyclic  loads  with 
maximum  nominal  stress  values  of  approximately  60  and  90  percent  of  the 
material  yield  strength  and  load  ratios  of  -1.0  and  0.1  with  a  frequency 
of  1.0  Hz.  The  numerical  simulations  of  specimens  consisting  of  Inconel 
71K  at  120<0^F~~("t  ime  -dependent  behavior)  consisted  of  cyclic  loads  where 
the  maximum  nominal  stress  was  .90  yield  strength,  the  applied  load 
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ratio  was  0.1,  and  the  cyclic  frequencies  were  .01  and  1.0  Hz. 

The  formation  of  a  plastic  wake  and  the  affects  of  plasticity 
induced  closure  were  observed  in  these  specimens.  Differences  in  the 
behavior  of  cracks  grown  by  cyclic  loading  and  those  exhibiting  no 
fatigue  growth  were  thus  investigated.  In  general,  it  was  found  that 


strain- related  characteristics  are  highly  dependent  upon  previous 


loading  history.  Stress-related  characteristics,  however,  were  found  to 


be  relatively  insensitive. 
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I.  INTRODUCTION 


Background 

The  United  States  Air  Force  places  great  emphasis  on  the  principles 
of  fracture  mechanics  in  determining  both  aircraft  construction  and 
maintenance  requirements  [1-3].  The  Engine  Structural  Integrity  Program 
(ENSIP)  provides  standards  for  structural  performance,  design 
development,  and  verification  requirements  for  turbine  engines.  This 
program  assumes  the  existence  of  initial  material  defects  which  may  be 
found  by  inspection  or  believed  to  exist  based  on  statistical  analyses. 
The  subsequent  prediction  of  material  life  considering  these  initial 
defects  is  an  important  aspect  of  this  program.  This  emphasis  on 
fracture  mechanics  also  incorporates  the  use  of  Non-Destructive 
Examinations  (NDEs)  of  aircraft  components  on  a  periodic  basis  to  detect 
flaws,  and  Linear  Elastic  Fracture  Mechanics  (LEFM)  principles  to 
predict  a  remaining  useful  component  life.  Using  these  methods,  if  a 
given  component  is  examined  and  it  is  determined  that  no  flaws  exist 
which  can  grow  to  a  critical  size  prior  to  the  next  periodic  examination 
or  scheduled  maintenance,  then  it  is  returned  to  service.  This  practice 
thereby  allows  for  an  efficient  use  of  aircraft  components  which  contain 
fatigue  damage  through  Retirement  for  Cause  (RFC)  procedures. 

One  possible  obstacle  to  these  approaches,  however,  is  the 
existence  of  short  cracks  and  their  apparent  non- conservative  growth 
patterns  [4-20] .  These  flaws  have  exhibited  growth  at  stress  levels  and 
sizes  when  LEFM  analyses  predict  that  no  crack  growth  would  occur. 
Furthermore,  the  growth  rates  of  these  short  cracks  are  significantly 
higher  than  that  predicted  using  LEFM  principles. 
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The  definition  of  a  short  crack  is  dependent  on  material,  geometry, 
crack  shape  and  other  possible  factors  [15].  It  is  therefore  not  the 
same  for  all  conditions.  In  general,  then,  a  crack  may  be  considered 
short  when  its  growth  rate  cannot  be  predicted  by  LEFM  principles. 

These  crack  sizes,  however,  are  less  than  the  crack  detection  limits  of 
current  NDE  and  inspection  methods  [15].  The  immediate  implication  is 
that  current  procedures  based  on  LEFM  principles  remain  valid.  So  long 
as  the  flaws  identified  by  these  procedures  remain  in  the  'long  crack' 
regime,  LEFM  analyses  will  continue  to  predict  fatigue  crack  growth 
accurately.  As  non -destructive  examination  and  inspection  methods 
improve,  however,  the  crack  detection  limit  will  likely  be  reduced  to 
the  point  where  short  cracks  can  be  identified,  and  the  non-conserva¬ 
tive  growth  of  these  flaws  will  become  a  major  factor  in  these 
procedures.  An  analysis  of  detectable  flaws  using  linear  elastic 
fracture  mechanics  principles  could  then  significantly  underestimate 
fatigue  crack  growth  and,  in  turn,  overestimate  an  aircraft  component’s 
serviceable  life. 

The  overall  consequence  of  this  situation  is  two-fold.  First,  the 
use  of  multi-million  dollar  aircraft  and  the  increasing  cost  of 
components  indicate  a  need  for  a  more  accurate  prediction  of  serviceable 
life.  A  non- conservative  estimate  based  on  LEFM  principles  could 
possibly  result  in  the  loss  of  an  aircraft  at  a  substantial  financial 
penalty  for  the  United  States  Air  Force.  Second,  and  most  importantly, 
is  the  safety  factor  involved  and  possible  loss  of  human  life  due  to  the 
overestimation  of  component  life. 

A  need  therefore  exists  to  understand  short  crack  behavior  and  the 
possible  factors  affecting  the  fatigue  growth  of  such  flaws.  Based  on 
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previous  studies  [21-23],  the  anomolous  behavior  of  short  cracks  may  be 
attributed  to  several  factors.  One  involves  microstructural 
considerations  such  as  grain  size  and  the  resulting  violation  of  the 
continuum  assumption  of  fracture  mechanics.  A  second  involves  the 
violation  of  LEFM  assumptions,  i.e.,  the  plastic  zone  sizes  of  short 
cracks  are  relatively  large  when  compared  to  their  lengths. 

It  has  also  been  suggested  that  the  short  crack  fatigue  response 
may  be  due  to  differences  in  the  crack  closure  response  of  long  and 
short  cracks  [22-34],  Crack  closure  may  be  envisioned  as  the  premature 
contact  of  crack  surfaces  during  the  unloading  portion  of  a  fatigue 
cycle  and  is  manifested  in  three  mechanisms:  asperity,  oxide,  and 
plasticity  induced  closure.  Both  asperity  induced  and  oxide  induced 
closure  involve  an  interference  between  the  crack  surfaces  which  prevent 
them  from  completely  closing  when  the  applied  load  is  reduced  to  zero. 

In  the  first  case,  this  interference  is  due  to  surface  roughness  or 
asperities,  whereas  in  Che  second  case,  it  is  due  to  oxides  forming  on 
the  crack  surfaces.  The  third  form,  plasticity  induced  closure,  is 
concerned  with  the  formation  of  plastically  deformed  material  at  and 
behind  the  crack  tip.  As  the  cracked  material  is  subjected  to  positive 
tensile  loading,  plastic  deformation  occurs  in  the  region  surrounding 
the  crack  tip  as  the  material  yield  strength  is  exceeded.  During 
subsequent  unloading,  the  surrounding  elastic  material  compresses  the 
plastic  region  causing  the  crack  to  close  prematurely.  It  is  this 
closure  mechanism  and  its  effect  on  short  crack  behavior  which  are  the 
areas  on  which  this  research  will  concentrate. 

With  a  thorough  knowledge  of  the  short  crack  phenomenon,  the  entire 
spectrum  of  crack  growth,  from  initiation  to  failure,  will  be  better 
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understood.  The  present  research  is  thus  aimed  at  developing  a  better 
understanding  of  short  crack  fatigue  growth  which,  in  turn,  will  lead  to 
more  accurate  predictions  of  component  serviceable  life  throughout  all 
ranges  of  crack  sizes.  A  more  efficient  and  cost  effective  use  of 
aircraft  components  and  Air  Force  resources  is  ultimately  envisioned. 

Literature  Review 

A  vast  amount  of  research  has  been  accomplished  and  literature 
published  in  the  attempt  to  understand  and  explain  short  crack  behavior. 
Recently,  a  number  of  works  [17,35,36]  have  been  published  discussing 
the  phenomenon  of  crack  closure  and  its  relation  to  short  crack  growth. 
Elber  [31]  found  that  closure  can  occur  during  the  tensile  portion  of 
the  fatigue  cycle  due  to  the  presence  of  residual  displacements  in  the 
crack  wake.  This  plasticity  induced  closure  results  from  the  constraint 
existing  around  the  wake  of  the  crack  tip  by  material  which  has  been 
permanently  deformed  within  the  plastic  zone  which  then  leads  to 
interference  between  the  opposing  crack  surfaces. 

Some  efforts  attempted  to  determine  when  and  why  LEFM  analyses 
become  non-conservative  [37-39]  based  upon  observations  of  short  crack 
growth  data.  Others  used  plasticity-adjusted  LEFM  or  elastic -plastic 
fracture  mechanics  to  describe  short  crack  behavior  [5,40,41],  Another 
group  considered  analytical  models  which  incorporated  the  effects  of 
crack  closure  and  residual  crack  tip  plasticity  [42-48], 

In  this  first  group  of  studies,  Allen  and  Sinclair  [37]  observed 
that  the  limitations  of  linear  elastic  fracture  mechanics  in  describing 
short  crack  behavior  may  be  due  to  an  inadequacy  in  describing  the 
stress  field  as  the  crack  length  approaches  zero.  A  nominal  plastic 


zone  size  was  calculated  assuming  that  the  high  near-tip  stresses  were 

redistributed  over  a  zone  through  yielding  and  that  this  zone  carries 

the  same  load  that  it  would  have  in  the  absence  of  yielding.  A  value  of 

K  ff,  that  stress  intensity  factor  which  produces  the  nominal  plastic 

zone  size,  was  then  determined.  With  this  value  of  K  conventional 

eff 

LEFM  was  used  to  satisfactorily  predict  threshold  stress  intensity 
factors  (those  values  of  Kj  below  which  a  crack  does  not  grow)  for  short 
flaws  in  center  cracked,  infinite  width  specimens  of  mild  and  4340 
steel . 

Hobson  [38]  proposed  a  crack  growth  equation  for  short  cracks  which 
could  be  used  in  conjunction  with  LEFM  to  describe  crack  growth  from 
crack  initiation  to  failure  in  an  aluminum  alloy.  Three  zones  of  crack 
growth  behavior  were  identified:  a  short  crack  zone  where  the  crack 
grows  to  a  length  corresponding  to  the  threshold  stress  intensity,  an 
interactive  zone  where  characteristics  of  both  long  and  short  crack 
growth  are  observed,  and  a  long  crack  zone  where  LEFM  principles  apply. 
'Best  fit'  analyses  were  then  accomplished  to  correlate  crack  growth 
data  with  the  proposed  equations. 

Taylor  and  Knott  [39],  in  a  correlative  analysis,  identified  three 
similar  regions  of  crack  growth  behavior  based  on  the  concept  of  a 
critical  length.  Below  this  value,  cracks  grow  faster  than  that 
predicted  by  LEFM.  Two  values  of  crack  length  (1^  and  lj)  which  are 
points  of  deviation  from  constant  stress  and  constant  stress  intensity 
behavior,  respectively,  were  determined.  At  crack  lengths  of  1^  and 
below,  it  was  found  that  the  initial  crack  length  had  no  effect  on 
fatigue  strength.  Below  values  of  1^,  the  LEFM  or  stress  intensity 
approach  becomes  non-conservative  in  predicting  fatigue  failure.  For  a 
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wide  range  of  materials,  a  correlation  between  the  values  of  1^  and  size 
of  the  material  microstructure  was  found. 

Whereas  these  first  efforts  attempted  to  correlate  short  crack 
behavior  with  LEFM,  the  second  group  used  plasticity-adjusted  LEFM  or 
elastic-plastic  fracture  mechanics  principles  to  describe  the  short 
crack  phenomenon. 

Trantina  and  de  Lorenzi  (41]  used  a  two-dimensional  finite  element 
model  to  analyze  short  crack  behavior  in  a  tension  loaded  double -edge 
cracked  specimen.  Plastic  flow  was  described  using  the  Ramberg- Osgood 
stress-strain  law  with  a  hardening  exponent  of  10  to  evaluate  the 
effects  of  elastic-plastic  material  behavior  on  the  stress  intensity 
factor.  J- integral  calculations  by  contour  integrals  and  virtual  crack 
extension  were  then  used  to  determine  an  effective  stress  intensity 
factor  value.  Cyclic  stress  versus  crack  length  data  showed  that 
deviations  from  LEFM  began  when  applied  stresses  were  about  0.7  of  the 
material  yield  strength  and  the  plastic  zone  size  was  approximately  0.1 
of  the  crack  length. 

El  Haddad  [5,40]  modified  both  elastic  and  elastic -plastic  fracture 
mechanic  solutions  to  correlate  short  crack  behavior.  By  introducing  an 
effective  length,  1^,  into  relations  for  the  stress  intensity  factor  and 
J- integral  and  by  expressing  the  stress  intensity  factor  in  terms  of 
strain  and  length  rather  than  stress  and  length,  El  Haddad  was  able  to 
adequately  correlate  short  crack  behavior.  Growth  rate  data  for  short 
flaws  in  notched  specimens  of  several  materials  correlated  well  with 
long  crack  data  under  this  approach. 

The  number  of  analytical  investigations  aimed  at  modeling 
plasticity  induced  closure  have  been  limited.  In  the  third  group  of 


6 


studies,  Newman  [42,44]  used  an  analytical  model  of  crack  closure  to 
study  crack  growth  behavior  and  closure  in  small  cracks.  The  closure 
model,  based  on  the  Dugdale  model  [43],  was  modified  to  leave 
plastically  deformed  material  along  the  crack  surface.  A  center 
cracked,  tension  specimen  and  symmetric  cracks  growing  from  circular 
holes  were  then  analyzed.  Crack  opening  stresses  calculated  using  this 
model  were  then  used  to  determine  effective  plastic  zone  corrected 
stress  intensity  factors  and,  in  turn,  applicable  crack  growth  rates. 
Newman  [44]  used  constant  strain  triangles  in  a  two-dimensional  finite 
element  analysis  which  accounted  for  non-linear  material  behavior  as 
well  as  changing  boundary  conditions  under  cyclic  loading.  To  satisfy 
these  boundary  conditions,  springs  were  attached  to  each  boundary  node, 
one  to  satisfy  conditions  in  the  x,  or  horizontal,  direction  and  one  in 
the  y,  or  vertical,  direction.  For  free  nodes,  the  spring  stiffness  was 
set  equal  to  zero,  whereas  for  fixed  or  closed  nodes,  it  was  set  to  some 
large  value.  The  particular  spring  stiffness,  depending  on  the 
appropriate  boundary  condition,  was  then  used  to  modify  the  elastic 
stiffness  matrix  as  required.  Assuming  elastic -perfectly  plastic 
material  under  constant  amplitude  or  two-level  block  loading,  crack 
growth  rates  were  found  consistent  with  experimental  data. 

Other  finite  element  analyses  have  been  used  to  analyze  crack 
closure  but  have  been  complicated  and  computer  intensive.  Ogura,  Oh j  i , 
and  Ohkibu  [45]  and  Nakagaki  and  Atluri  [46]  utilized  higher  order 
finite  element  models  to  study  crack  closure  in  various  specimens. 

Based  upon  the  conclusions  found  in  the  preceeding  literature,  as 
well  as  those  in  other  works  [47-62],  several  items  are  indicated. 

First,  the  study  of  short  cracks  requires  an  elastic-plastic  analysis. 
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LEFM  analyses  do  not  accurately  predict  short  crack  behavior.  Second, 
numerous  analytical  models  have  been  developed  which  can  satisfactorily 
predict  short  crack  behavior  or  closure  in  fatigued  specimens.  And 
third,  numerical  models,  primarily  finite  elements,  have  been  employed 
to  either  predict  short  crack  behavior  or  crack  closure  effects.  None, 
however,  have  been  found  which  combine  the  two  phenomena. 

As  a  result,  this  research  concentrates  on  a  finite  element 
analysis  of  plasticity  induced  closure  and  its  effect  on  short  crack 
fatigue  growth.  Initially  short  cracks  growing  from  the  edge  of  a 
finite-width  specimen  under  constant  amplitude  fatigue  cycling  at 
various  stress  levels  and  stress  ratios  are  investigated.  Crack  closure 
is  traced  using  a  two-dimensional  elastic -plastic  finite  element 
analysis.  Results  are  interpreted  in  terms  of  the  influence  of 
plasticity  induced  closure  on  the  growth  of  small  fatigue  cracks. 
Observations  as  to  the  sensitivity  of  closure  to  the  measurement 
location  along  the  crack  length  are  also  discussed. 

Approach 

As  previously  stated,  this  research  concentrates  on  developing  a 
better  understanding  of  plasticity  induced  closure  and  its  effect  on 
short  crack  behavior  in  a  cyclically  loaded  specimen.  The  subject 
specimen  studied  here  is  composed  of  Titanium  6Al-2Sn-4Zr-6Mo  (TI-6246) 
at  room  temperature  and  considered  to  be  a  homogeneous,  isotropic 
material  that  deforms  and/or  strains  under  cyclic  loading.  Starting 
with  an  initial  single-edge  crack  (a^  -  .001  inch),  the  material  is 
subjected  to  various  load  ratios  and  amplitudes  in  order  to  investigate 
the  short  crack  phenomenon.  A  similar  specimen,  composed  of  Inconel  718 
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at  1200°F  is  also  studied  to  incorporate  viscoplastic  effects.  Two 
separate  load  frequencies  are  applied  to  isolate  this  phenomenon. 

Visco  II,  an  in-house  computer  program,  was  modified  to  incorporate 
the  effects  of  crack  closure.  This  finite  element  code  employs  constant 
strain  triangular  elements.  It  has  the  ability  to  apply  saw-tooth 
shaped  cyclic  loads  and  to  simulate  crack  growth  through  the  release  of 
finite  element  nodes.  In  addition,  the  Bodner- Partom  constitutive 
relations  [63-65]  are  employed  to  model  material  behavior  during  load 
increments.  This  elastic-viscoplastic  model  is  primarily  time- 
dependent  ,  but  can  easily  be  made  rate  independent  making  it  a  versatile 
tool  in  this  research  effort.  An  explicit  Euler  extrapolation  scheme  is 
used  to  Integrate  these  relations  with  respect  to  time  and  is  further 
incorporated  in  Visco  II  using  the  residual  force  technique. 

A  finite  element  mesh  was  developed  to  model  a  single-edge  cracked 

specimen  (a.  -  .001  inch,  W  -  .  10  inch)  and  was  investigated  to  compare 

the  effects  of  various  cyclic  loading  profiles.  For  the  TI-6246 

specimen,  external  load  was  input  as  a  saw-tooth  s-ress-tiine  pattern  of 

constant  amplitude  with  load  ratios  (ratio  of  minimum  to  maximum  load) 

of  -1.0  and  0.1.  Maximum  load  amplitudes  were  applied  to  obtain  nominal 

stresses  of  approximately  60  and  90  percent  of  the  material  yield 

strength.  These  resulted  in  initial  stress  intensity  factors  (K^)  of 

1/2 

6.69  and  10.04  KSI-IN  respectively.  Cyclic  load  frequency  in  chest- 

cases  was  maintained  at  1.0  Hz.  In  the  study  of  Inconel  718,  a  load 

ratio  of  0.1  and  a  nominal  maximum  stress  value  of  .90  yield  strength 

were  maintained.  The  initial  stress  intensity  factor  for  this  specimen 
1/2 

was  7.354  KSI-IN  However,  the  cyclic  frequency  took  on  values  of 

1 . 0  and  .01  Hz . 
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II.  Theory 


In  this  research  effort,  a  single-edge  cracked  specimen  (a^  -  .001 
inch)  is  subjected  to  cyclic  loading  in  order  to  study  the  effects  of 
crack  closure  and  the  formation  of  a  plastic  wake  as  a  short  crack 
propagates.  Using  a  finite  element  analysis,  specimens  composed  of 
materials  displaying  entirely  different  characteristics  (TI-6246  at  room, 
temperature  and  Inconel  718  at  1200°F)  are  investigated.  Inconel  718 
exhibits  non-linear,  viscoplastic  (time-dependent)  material  behavior  at 
high  temperatures,  whereas  TI-6246  at  room  temperature  is  considered 
non-linear,  elastic -plastic  ( time  -  independent) .  Due  to  the  variability 
of  the  materials  involved,  a  viscoplastic  flow  law  is  required  which  can 
easily  be  modified  to  account  for  time  independent  relations.  In 
addition,  a  finite  element  code  is  needed  which  incorporates  cyclic 
loading,  fatigue  crack  growth,  and  crack  closure.  This  section 
discusses  the  applicable  relations  chosen  to  model  plastic  flow  as  well 
as  finite  element  considerations  required  for  such  an  analysis. 

Bodner- Partom  Viscoplastic  Flow  Law 

The  Bodner-Partom  Viscoplastic  Flow  Law  (Bodner  model)  satisfies 
the  above  requirements  for  both  time -dependent  and  time -independent 
behavior  and  has  been  used  successfully  in  the  past  to  study 
viscoplasticity  in  a  variety  of  analyses  [66-71]  including  cyclically 
loaded  specimens.  This  law,  based  on  dislocation  dynamics,  suggests  a 
continuous  flow  relation  between  stress  and  viscoplastic  strain 
commencing  at  the  onset  of  loading.  As  used  here,  this  law  assumes  that 
the  materials  involved  are  isotropic  and  exhibit  no  kinematic  hardening. 
The  following  discussion  is  intended  to  provide  a  basic  understanding  of 
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the  Bodner  model  and  present  key  items  pertinent  to  this  research.  More 
detailed  information  may  be  obtained  in  works  directly  related  to  the 
model  [63-65] . 

Formulation  of  the  constituitive  equations  begins  by  recognizing 
that  for  small  strains,  total  strain  can  be  separated  into  elastic,  as 
well  as  viscoplastic,  strains  such  that 


€  .  . 
ij 


(2.1) 


Taking  the  time  derivative  of  this  equation,  the  total  strain  rate  is 
then  expressed  as 


+  € 


ij 


(2.2) 


•  0 

where  the  elastic  strain  rate,  €. .  ,  is  related  to  the  stress  rate 

ij 

through  the  time  derivative  of  Hooke's  law. 

The  determination  of  the  viscoplastic  strain  rate,  however,  begins 
by  squaring  the  rate  form  of  the  Prandtl -Reuss  plastic  flow  equation, 


€.  . 

IJ 


P 


(2.3) 


to  obtain 


k.  .  P€  ,  . P  -  x2s.  S  . 

ij  ij  ij  ij 


(2. A) 


where  €.jP  are  components  of  the  deviatoric  viscoplastic  strain  rate 
tensor,  S„  are  components  of  the  deviatoric  stress  tensor,  and  X  is  a 
proportionality  constant.  Then  by  substituting  the  following 
relations , 


J0  -  (S.  .S. . )/2 

2  ij  ij" 


(2.5) 
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(2.6) 


and  D2P  ‘  <  VV>/2 

one  obtains 


(2.7) 


Here,  is  defined  as  the  second  invariant  of  the  plastic  strain  rate 

and  is  the  second  invariant  of  the  deviatoric  stress  tensor. 

Based  on  experimental  work,  Bodner  and  Partom  expressed  D2P  as 


D2P  -  D02exp[-(Z2/3J2)n((n+l)/n)] 


(2.8) 


In  this  relation,  Dg  is  the  limiting  value  of  the  plastic  strain  rate  in 
shear,  Z  is  a  measure  of  material  hardness,  and  n  a  rate  sensitivity 
parameter.  Further,  the  evolution  equation  for  Z  is  assumed  to  consist 
of  a  hardening  term  due  to  plastic  work,  Wp,  and  can  be  described  as 

Z  -  (dZ/fcU  )(dW  /dt)  (2.9) 

P  P 

Here,  (3Z/&W  )  is  assumed  to  be  of  the  form 
P 

OZ/fcWp)  -  m(Zx-Z)  (2.10) 

where  m  is  a  parameter  controlling  the  rate  of  material  work  hardening 

and  Z.  is  the  maximum  value  of  Z  attainable.  W  is  the  relative  amount 
1  P 

of  plastic  work  done  from  some  initial  state  and  is  expressed  in  the 
form 

“p  '/SijlijPdt  '  2<D2Pj2)1/2  (2  U) 

Most  materials  exhibit  some  thermal  recovery  of  hardening  or 
relaxation  of  accumulated  plastic  work  at  high  temperatures.  To  model 
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this  particular  behavior,  the  change  in  hardness  can  be  modified  by 
adding  the  term 

Z  -  -A((Z-Z0)/Z. )rZ.  (2.12) 

rec  2'  11 

In  this  relation,  Z2  is  the  value  of  Z  in  a  completely  non-work  hardened 
condition  indicating  the  minimum  value  of  material  hardness,  whereas  A 
and  r  are  the  hardening  recovery  coefficient  and  exponent  respectively. 
Then,  by  integrating  the  resulting  equation  for  the  change  in  hardness 
using  equation  (2.10),  material  hardness  takes  on  the  form 

Z  -  Zx  +  (ZQ-Z1)exp(-mWp/Z0)  -  AUZ-Z^/Z^1^  (2.13) 

where  Zq  is  defined  as  the  initial  value  of  hardness  prior  to  any 
plastic  deformation. 

For  strictly  time  independent  behavior,  the  thermal  recovery  term 

in  equation  (2.13)  is  negligible  and  can  therefore  be  disregarded  in  the 

evaluation  of  material  hardness.  This  can  be  accomplished  by  setting 

the  values  of  A,  r,  and  Z0  to  zero  in  the  Z  term  without  affecting 

2  rec  ° 

the  overall  solution  of  the  Bodner  equations. 

The  specific  procedure  used  to  solve  the  Bodner  equations  is  an 
Euler  extrapolation  scheme  which  integrates  the  equations  with  respect 
to  time.  This  is  accomplished  for  each  element  in  the  finite  element 
code  as  follows: 


Z1  -  Z1  +  (Z0-Z1)exp[-»Wp1'1/Z0] 

(2.14) 

(D2V  -  D02exp[-((Z1)2/3J2i‘1)n((n+l)/n)] 

(2.15) 

(^ijP)1  ~  KD2P)i/J21'1]1/2(Sij)1*1 

(2.16) 
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(2.17) 


up‘  -  “p1'1  *  lsij|1'1|d£ijP|1  -  ^/otVl^z^z1)] 

where  the  superscript  i  refers  to  the  current  timestep. 


(2.18) 

(2.19) 


Material  Constants 

The  constants  used  in  the  Bodner  model  are  both  material  and 
temperature  dependent  and  must  therefore  be  evaluated  for  each  material 
at  the  specific  temperature  of  interest.  The  Bodner  constants  for 
Inconel  718  at  1200°F  have  previously  been  determined  by  Beaman  [72]  and 
are  presen“ed  in  Table  2.1.  These  values  were  subsequently  verified  by 
Mercer  [73]  and  have  been  shown  to  model  the  particular  material 
behavior  satisfactorily. 

Material  constants  for  TI-6246  at  room  temperature  have  not  been 
previously  evaluated  for  use  with  the  Bodner  model  and  must  therefore  be 
determined  from  available  experimental  data.  Following  the  procedures 
developed  by  Beaman  [72]  and  Stouffer  [74]  in  conjunction  with 
experimental  stress  -  strain  data  provided  by  the  WRDC  Materials 
Laboratory,  values  for  these  Bodner  constants  were  determined  and  are 
listed  in  Table  2.2.  An  overview  of  the  steps  involved  in  this  approach 
is  included  in  Appendix  A.  Note  that  the  parameters  affecting  thermal 
recovery  (A,  r,  and  2^)  have  been  given  the  value  of  zero.  TI-6246  at 
room  temperature  exhibits  time  independent  behavior  and  the  thermal 
recovery  term  in  equation  (2.13)  is  negligible. 

The  values  of  the  Bodner  constants  determined  for  TI-6246  at  room 
temperature  were  subsequently  verified  by  analytically  loading  a 
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Table  2.1 


Bodner  Coefficients  for  Inconel  718  at 


Parameter 


Description 


Elastic  Modulus 


23 

(16 


Strain  Rate  Exponent 


Limiting  Value  of 
Strain  Rate 


Limiting  Value  of 
Hardness 


Maximum  Value  of 
Hardness 


Minimum  Value  of 
Hardness 


Hardening  Rate  Exponent 


Hardening  Recovery  Coefficient  1 


Hardening  Recovery  Exponent 


(1  Kbar  -  100  MPa  -  14.504  KSI) 


12006F 


Value 


.568  x  10^  KSI 
.250  x  10A  MPa) 


3.0 


m6  ' 1 

10  sec 


235.3  KSI 
(1622  MPa) 

260.3  KSI 
(1795  MPa) 

104.1  KSI 
(718  MPa) 

2.875  KSI 'l 
(.417  MPa  ) 

.5  x  10  ^  sec 

7.00 
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Table  2.2 


Bodner  Coefficients  for  TI-6246  (Room  Temperature) 


Parameter 


Description 


Value 


Elastic  Modulus 


16.275  x  10"  KSI 
(11.221  x  10  MPa) 


n 


Strain  Rate  Exponent 


2.0 


Limiting  Value  of 
Strain  Rate 


lf,4  -1 
10  sec 


Limiting  Value  of 
Hardness 


351.72  KSI 
(2425.0  MPa) 


Maximum  Value  of 
Hardness 


396.87  KSI 
(2736.3  MPa) 


Minimum  Value  of 
Hardness 


0.0 


m 


Hardening  Rate  Exponent  .069  KSI. 

/  Ai  wrt  "  ^ 


-1 


(.01  MPa  ) 


Hardening  Recovery  Coefficient  0.0 


Hardening  Recovery  Exponent 


0.0 


(1  Kbar  -  100  MPa  -  14.504  KSI) 
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uniaxial,  two-element  finite  element  model  under  a  similar  stress  range 
and  approximate  strain  rate  of  .005  in/min.  Applied  loads  ranged  from  0 
to  176  KSI  to  coincide  with  the  experimental  data.  The  Visco  code, 
however,  is  not  governed  by  strain  rate  but  instead  uses  a  load  versus 
time  input  format.  The  strain  rate  of  .005  in/min  was  approximated  by 
varying  the  rate  at  which  the  load  was  applied  until  a  strain  rate 
approximating  .005  in/min  was  obtained.  The  experimental  stress-strain 
data  and  the  results  from  the  two-element  finite  element  model  are 
presented  in  Figure  2.1.  Since  Visco  is  load  or  stress  driven,  strain 
rates  cannot  be  maintained  during  cyclic  loading.  For  this  reason, 
further  analysis  of  TI-6246  material  constants  could  not  be  performed. 
The  results  of  the  two-element  finite  element  model,  therefore,  were 
deemed  sufficient  to  verify  the  Bodner  constants  for  this  research. 

Finite  Element  Considerations 

Visco  II,  the  computer  program  used  in  this  research,  is  a 
two-dimensional,  plane  stress/strain,  finite  element  code  which  uses 
constant  strain  triangular  elements.  The  overall  solution  technique  is 
the  residual  force  method  and  key  features  incorporated  within  this 
program  are  a  Gauss-Siedel  iterative  procedure  coupled  with  a  crack 
closure  algorithm  to  solve  for  equilibrium,  the  Euler  extrapolation 
scheme  previously  discussed  to  model  elastic-viscoplastic  material 
behavior,  a  data  storage/restart  capability  to  handle  extensive  cyclic 
loading,  and  a  procedure  to  simulate  crack  growth.  It  should  be  pointed 
out  that  Visco  II  is  an  extension  of  the  Visco  program  developed  by 
Hinnerichs  [71,75],  modified  to  meet  the  needs  of  this  particular 
research  effort. 
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250. 


Figure  2.1 

Stress-Strain  Data  -  TI-6246 
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.016  .024  .032  .040  .048  .056 

y-STRRIN  UN/IN) 


Residual  Force  Method.  The  basic  concept  in  the  residual  force 
method  is  that  the  elastic  stiffness  matrix  remains  constant  throughout 
an  analysis  and  plasticity  effects  are  incorporated  into  the  solution 
through  the  inclusion  of  a  plastic  load  vector.  The  resulting 
equilibrium  equations  then  take  on  the  form: 

[K] (u) 1  -  { P) 1  +  (Q } 1  (2.20) 

where  [ K ]  is  the  constant  elastic  stiffness  matrix,  |u)^  is  the  nodal 
displacement  vector,  |P)^  is  a  vector  of  externally  applied  loads,  and 
IQ)1  is  the  vector  of  internal  loads  developed  due  to  the  accumulation 
of  plastic  deformation.  The  superscript  i  in  this  relation  refers  to 
the  current  timestep.  The  algorithm  employed  in  Visco  II,  then,  takes 
on  the  following  form: 

1.  Compute  the  current  time, 

t1  -  t1'1  +  dt1  (2.21) 


(2.22) 

(2.23) 

(2.24) 


2.  Calculate  the  plastic  strain  rate 
using  the  Bodner  model/Euler 
extrapolation  scheme  (see  Eq  2.16), 


<€.  .V  -  [(D2P)i/J2i'1]1/2(S.  .  I1'1 


3.  Compute  the  plastic  strain  increments, 


<Vl  -  1 V1-'1 


4.  Calculate  the  plastic  load  vector  (elementwise), 

IQ! 


B]T[D](CijP)idvol 


19 


where 


l€ijPli  “  'V*1'1  +  ,d€ijP)1 


5.  Compute  the  current  external  loads, 


{ P} 1  -  (P)1'1  +  { P> idti 


Calculate  nodal  displacements  using 
the  Gauss-Siedel  technique, 

N 


(u)1  -  IP)1  +  IQ)1) 

J=1 


If  closure  has  occurred,  determine  the 
resultant  nodal  forces  and  modify  the 
external  load  vector  accordingly, 


(P}1  -  IP)1  +  (Pclo)i 


and  re -solve  for  nodal  displacements 


8.  Calculate  total  strains  for  each 
individual  element. 


(€ij)1  -  ( B]  ( u) 1 


9.  Compute  resultant  stresses, 


lo\j ) 1  -  [DKle.j)1  -  tc^V) 


10.  Modify  current  timestep  size 
and  continue  procedure  until 
analysis  is  completed. 


(2.25) 


(2.26) 


(2.27) 


(2.28) 


(2.29) 


(2.30) 
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With  this  procedure,  time  is  incremented  directly  and  the  resultant 
displacements,  stresses  and  strains  are  subsequently  evaluated.  The 
particular  timestep,  dt,  is  determined  subject  to  stress  and  strain 
tolerances  and  an  algorithm  relating  these  constraints  to  changes  in  the 
current  timestep  size  [75-77].  This  approach  thereby  maximizes  the 
allowable  timestep  and  thus  minimizes  computational  requirements.  The 
stress  and  strain  parameters,  Pg.  and  P£ ,  used  for  developing  timestep 
sizes,  are  functions  of  effective  stress  and  effective  strain  and  are 
defined  as  follows  [75-77]: 


(2.31) 


P£  -  <0  £t0l 


(2.32) 


Pa  and  P£  are  evaluated  for  each  element  during  each  timestep  and  the 
parameter  P  is  set  equal  to  the  largest  value.  The  timestep  is  then 
determined  according  to  a  generalized  set  of  relationships: 

dt*  -  c^dt*  *  if  <  P  <  Cj  (2.33) 

where  c^  is  a  proportionality  constant,  and  the  values  of  c ^  and  c^  form 
intervals  for  the  value  of  P.  By  controlling  the  size  of  subsequent 
timesteps  in  this  manner,  the  levels  of  accuracy  inherent  in  the 
equilibrium  and  compatability  relations  are  also  determined.  The 
evaluation  of  these  constants  is  discussed  in  a  following  section 
related  to  finite  element  modeling. 


Gauss-Siedel 


The  equilibrium  equations  are  solved  in 


Visco  II  using  a  Gauss-Siedel  iterative  technique  [75],  Nodal 
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displacements  are  determined  individually  and  assumed  to  be  a  function 
only  of  those  nodes  directly  connected  to  the  one  in  question.  At  the 
completion  of  each  iteration,  the  change  in  nodal  displacements  is  used 
to  calculate  a  tolerance  limit  which  is  then  used  to  check  for 
convergence.  This  is  dependent  upon  the  initialization  of  nodal 
displacements  which  is  normally  set  to  zero.  However,  displacements  from 
previous  solutions,  representing  a  more  accurate  'guess',  will  speed  the 
convergence  of  this  technique.  Additionally,  an  overrelaxation  factor 
is  employed  in  order  to  speed  the  process  [71]. 

If  negative  displacements  occur  along  the  crack  face  indicating 
crack  closure,  Visco  II  performs  an  iterative  closure  routine.  This 
procedure  first  resets  applicable  nodal  displacements  to  zero  (as 
negative  values  are  physically  impossible),  and  then  calculates  a 
resultant  nodal  force  required  to  insure  a  zero  displacement  using  the 
following  relation: 


(P 


clo 


i 


T 


l<7 


.  .  )  dvol 

ij 


(2.34) 


{P  ,  l1  is  the  vector  of  external  forces  at  the  closed  nodes  and  is  a 
clo 

function  of  stresses  in  the  three  adjacent  elements.  These  forces  are 
then  added  to  the  external  load  vector  and  the  solution  for  equilibrium 
is  repeated.  This  procedure  thus  incorporates  the  changing  boundary 
conditions  involved  in  crack  closure  without  modifying  the  stiffness 
matrix  as  is  the  case  in  other  closure  routines  [42,47,50].  Since  nodal 
displacements  will  change  with  the  addition  of  these  closure  loads,  this 
crack  closure  algorithm  is  repeated  after  the  equilibrium  solution  until 
the  difference  in  closure  loads  on  subsequent  iterations  is  less  than  a 
specified  tolerance  value.  This  procedure  can  be  visualized  in  the 
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i>  NO  CLOSURE. 


M1  *  0 


ID  NEGATIVE  DISPLACEMENTS. 


ill)  RESET  DISPLACEMENTS  TO  ZERO. 

a  ■  a 


|y)  DETERMINE  CLOSURE  LOADS.  W  -  MM' 

M  -  B>We,P  -  M'\ 

•  2/(b]|cr,)i'dvol 

v)  RE-SOLVE  EQUIUBRIUM.  \u\'  -  [K]*'}|P{'  ♦  ♦  $Q('  ( 

Figure  2.2 

Finite  Element  Crack  Closure 
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Continue 


Figure  2 . 3 

Equilibrium/Closure  Solution  Algorithm 


diagram  in  Figure  2.2  and  flowchart  in  Figure  2.3.  Because  the 
procedure  is  only  performed  when  negative  displacements  occur  at  nodes 
along  the  crack  surface,  subsequent  crack  opening  is  easily  handled. 
Crack  opening  is  indicated  by  a  lack  of  negative  displacements  and  the 
crack  closure  algorithm  is,  thus,  not  performed.  Additionally,  by 
compensating  for  crack  closure  through  the  inclusion  of  an  internal  load 
vector,  the  global  stiffness  matrix  remains  constant  and  does  not 
require  any  modifications  during  the  above  routine. 

Crack  Growth .  In  Visco  II,  crack  growth  is  simulated  through  the 
use  of  a  node  release  algorithm  [71,75].  When  a  node  is  to  be  released, 
the  forces  at  that  node  required  to  maintain  a  zero  displacement  are 
evaluated  using  equation  (2.34).  This  is  equivalent  to  an  oppositely 
directed  force  and  is  then  applied  to  the  node  and  reduced  linearly  over 
time  (Figure  2.4).  The  boundary  condition  that  had  previously  been  in 
effect  is  conveniently  handled  in  the  Gauss-Siedel  technique  previously 
discussed.  The  stiffness  matrix  is  initially  developed  on  an  individual 
nodal  basis  in  which  only  adjacent  nodes  contribute  to  one's 
displacement.  When  a  node  is  fixed,  its  equilibrium  equation  is  skipped 
over  in  the  iterative  solution  procedure.  As  that  node  is  released, 
however,  its  equilibrium  equation  is  then  included.  By  handling  crack 
growth  and  changing  boundary  conditions  in  this  manner,  refactoring  of 
the  stiffness  matrix  is  not  required. 

Data  Storage /Restart  Capability.  Visco  II  further  includes  a 
routine  which  provides  for  the  ability  to  store  data  at  the  completion 
of  the  finite  element  analysis  and  subsequent  ability  to  restart  the 
program  using  this  data  as  initial  conditions.  In  this  manner,  it  is 
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possible  to  continue  an  extensive  analysis  indefinitely  without  being 
limited  by  computer  restraints.  Additionally,  this  capability  provides 
a  method  to  perform  a  finite  element  analysis  at  any  point  in  a  fatigue 
crack  growth  study,  the  only  requirement  being  that  initial  conditions 
are  known  or  available.  This  program  modification  is  outlined  in 
Appendix  B. 


— H  k—  node  release 

TIME 


Figure  2.U 

Node  Release  Method  (a)  Finite  Element  Cracktip 

(b)  Node  Unload  Method 
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Ill .  Finite  Element  Model 


In  the  previous  chapter  discussing  theory,  several  programming 
considerations  were  presented  which  could  be  used  to  solve  a  generalized 
elastic -viscoplastic  flow  problem  using  the  finite  element  method. 

These  were  incorporated  into  the  finite  element  code,  Visco  II,  and 
included  the  particular  viscoplastic  flow  law  (Bodner  model)  and 
associated  material  parameters,  solution  procedures  such  as  the  residual 
force  method  and  Gauss-Siedel  iterative  technique,  as  well  as  several 
other  considerations  relevant  to  this  research.  In  the  sections  to 
follow,  however,  more  specific  information  is  presented  which  is  a 
function  of  this  particular  research  effort.  Included  here  are  the 
investigation  and  selection  of  suitable  programming  variables,  valida¬ 
tion  of  program  subroutines,  and  the  determination  of  an  appropriate 
finite  element  mesh.  The  information  determined  in  these  areas  was 
based  on  a  series  of  preliminary  studies  conducted  as  an  integral  part 
of  this  research  effort  and  are  included  in  this  chapter  for 
completeness . 

Investigation  of  Programming  Variables. 

Within  the  residual  force  algorithm  employed  in  Visco  II,  the 
accuracy  of  the  equilibrium  and  compatibility  relations  is  controlled 
through  the  use  of  the  stress  and  strain  tolerances  coupled  with  a  set 
of  relations  to  modify  the  current  timestep  size  as  discussed  in  chapter 
2.  These  relations  are  of  the  form: 

dt1  -  Cjdt1  if  c2  <  P  <  c3  (2.33) 

where  P  is  the  greatest  value  of  P^.  and  P£  determined  by  equations 
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(2.31)  and  (2.32),  is  a  proportionality  constant,  and  c ^  and  form 
intervals  around  P.  In  order  to  determine  the  optimal  values  of  these 
parameters,  a  one-dimensional  case  was  investigated.  A  two-element 
model  was  subjected  to  a  uniaxial,  fully  reversed,  cyclic  load  in  which 
the  values  of  c^,  c^ ,  and  c^  were  varied  to  study  their  effects  on 
accuracy  and  computer  time  requirements,  the  results  of  which  are  shown 
in  Figures  3.1  and  3.2.  In  these  figures,  t(cps)  refers  to  the  computer 
time  required  to  complete  one  entire  load  cycle,  whereas  r  is  the  amount 
of  drift  or  ratcheting  exhibited  in  the  stress  -  strain  curve  over  the 
same  period.  The  values  of  'base  r'  and  'base  t'  are  those  initially 
obtained  using  the  relations  for  timestep  change  set  forth  by  Hinnerichs 
[71,75]  and  used  in  subsequent  studies  [67-70].  Based  upon  these 
results,  the  complete  set  of  relations  take  on  the  following  form: 


dt1  -  . 9dtX 
dt1  -  l.Odt1 
dt1  =  l.ldt1 
dt1  -  1 . 2dt* 


if  1.0  <  P 

if  .85  <  P  <  1.0 

if  .70  <  P  <  .85 

if  P  <  .70 


(3.1) 

(3.2) 

(3.3) 

(3.4) 


The  equilibrium  equations  are  solved  for  in  the  residual  force 
method  employed  here  using  the  Gauss- Siedel  iterative  technique.  These 
equation?;,  in  matrix  form, 

[K] lu)1  -  (P)1  +  (Q)1  (2.20) 

are  solved  for  the  nodal  displacements  during  each  timestep.  The  number 
of  iterations  required  to  solve  these  equations  is  controlled  by  a 
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Effect  of  Varying  in  Timestep  Algorithm 


Figure  3.2 


Effect  of  Varying  of  and  c3  in  Timestep  Algorithm 
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'force  unbalance'  term  which  represents  the  numerical  error  in  the 
equilibrium  solution.  At  the  completion  of  each  iteration,  the  force 
unbalance  is  determined  by  summing  the  forces  resulting  from  the 
incremental  displacements  calculated  during  that  iteration.  This 
procedure  continues  until  this  value  is  less  than  some  tolerance  value 
or  the  number  of  iterations  exceeds  a  predetermined  number  [71]. 

The  center  cracked  panel  employed  by  Henkel  [68]  was  subjected  to  a 

fully  reversed  cyclic  load  with  a  stress  intensity  factor  of  35  KSI- 
1  /2 

IN  '  to  investigate  this  concept.  In  this  study,  values  of  the  force 
unbalance  tolerance,  maximum  number  of  iterations  allowed,  and  the 
overrelaxation  factor  were  all  varied.  In  each  instance,  the  value  of 
the  force  unbalance  diminished  at  an  exponential  rate  until  additional 
iterations  did  not  affect  its  final  value.  Further,  changes  in  the 
force  unbalance  tolerance  and  overrelaxation  factor  did  not  substan¬ 
tially  improve  this  value.  This  study  showed  that  although  the  accuracy 
of  the  equilibrium  solution  could  not  be  significantly  improved  by 
varying  the  tolerance  values  within  the  Gauss-Siedel  technique,  the 
amount  of  computer  time  required  could  be  reduced  by  limiting  the  number 
of  iterations  allowed. 

The  stress  and  strain  tolerances  employed  in  the  residual  force 
method  were  further  investigated  by  subjecting  a  single-edge  cracked 
specimen  to  cyclic  loading  at  a  nominal  stress  value  of  approximately  60 
percent  material  yield  strength  and  a  load  ratio  of  0.1.  The  stress  and 
strain  tolerances,  a  ^  and  e  ^  respectively,  were  found  to  optimize 
the  tradeoff  between  accuracy  and  computational  requirements  at  a  value 
of  .10.  That  is,  the  amount  of  change  in  stress  or  strain  in  any 
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element  is  not  allowed  to  exceed  ten  percent  during  a  single  time 
increment.  If  this  amount  of  change  is  exceeded,  then  the  current 
timestep  is  reduced  and  the  solution  repeated.  At  tolerances  less  than 
ten  percent,  these  values  were  frequently  exceeded  which  resulted  in 
significantly  higher  computational  requirements.  Under  these 
circumstances  and  the  additional  requirement  to  analyze  the  specimens 
under  a  nominal  stress  of  .90  yield  strength,  these  particular  tolerance 
values  were  deemed  appropriate  for  the  analysis  conducted. 

Results  obtained  from  these  analyses  agreed  with  those  obtained  by 
Henkel  [68]  and  were  subsequently  verified  by  Mercer  [73]  in  an 
independent  study  using  an  implicit  solution  technique.  Representative 
strain  fields  under  similar  conditions  may  be  found  in  works  by  Wilson 
[70]  and  Henkel  [68]. 


Crack  Closure 

Within  the  Gauss-Siedel  technique  employed  in  Visco  II  is  the  crack 
closure  algorithm  previously  discussed  in  chapter  2.  When  closure 
occurs,  a  resultant  internal  force  using 


(P 


clo 


1  -L/*BiT(aij,dvo1 


(2.34) 


is  determined  and  applied  to  the  closed  node  by  modifying  the  external 
load  vector: 


IP)1  - 


IP)1  + 


{Pclo> 


(2.28) 


In  this  manner,  the  global  stiffness  matrix  does  not  require 
modification  and  thus  maintains  the  overall  solution  technique  employed 
by  the  residual  force  method. 
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This  closure  algorithm  was  applied  to  the  center  cracked  panel 

investigated  by  Henkel  [68]  using  a  stress  intensity  factor  of  45 

KSI-IN*^  and  load  ratio  of  -1.0.  In  Henkel's  work,  where  the 

displacements  at  closed  nodes  were  simply  reset  to  zero,  residual  force 

unbalances  on  the  order  of  2000-3000  were  evident  at  the  point  of  full 

negative  load.  Through  the  use  of  this  closure  algorithm,  the 

magnitudes  of  the  force  unbalance  were  reduced  to  less  than  100  at  the 

same  load  points.  By  thus  reducing  the  residual  force  unbalance  in  the 

Gauss-Siedel  equilibrium  solution,  the  overall  accuracy  was  improved. 

Further  verification  of  this  routine  was  accomplished  by  analyzing  a 

1/2 

center  cracked  panel  with  a  stress  intensity  factor  of  35  KSI-IN 
under  fully  reversed  loading.  The  results  of  this  analysis  were 
compared  to  those  obtained  by  Mercer  [73]  for  the  identical  specimen 
using  a  completely  different  finite  element  code.  Stresses  obtained 
from  these  two  programs  for  elements  ahead  of  the  crack  tip  are 
displayed  in  Figure  3.3 

Crack  Growth 

Crack  growth  is  simulated  using  the  node  release  method  previously 
discussed  and  has  been  successfully  used  in  several  past  finite  element 
analyses  [69,71,75],  Node  release  (crack  growth)  is  assumed  to  occur  at 
the  point  of  maximum  load  application  ,i.e.,  at  the  one-quarter  cycle 
point  for  R  -  -1.0  and  at  the  one-half  cycle  point  for  R  -  0.1.  This 
allowed  the  forces  to  be  reduced  at  the  released  node  as  the  total 
external  loads  were  also  being  reduced,  the  result  being  that  stresses 
near  the  new  crack  tip  did  not  increase  significantly  during  the  node 
release  process.  This  condition  may  have  resitted  had  node  release 
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Figure  3.3 

Crack  Closure  Algorithm  Results 
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occurred  during  an  increase  in  external  loads.  As  explained  in  chapter 
2,  when  a  node  is  released,  a  force  is  initially  calculated  and 
internally  applied  to  that  node  to  maintain  a  zero  displacement.  This 
force  is  then  reduced  over  a  specific  time  interval,  the  node  release 
time,  until  that  node  is  no  longer  restrained.  For  the  finite  element 
simulation  here,  the  node  release  time  was  arbitrarily  set  equal  to 
one-half  the  time  from  maximum  load  application  to  the  point  of  minimum 
positive  load.  Hence,  for  R  -  -1.0,  this  time  was  set  at  .125  seconds, 
whereas  the  release  time  was  .25  seconds  for  a  load  ratio  of  0.1.  In 
this  manner,  the  node  release  algorithm  was  completed  prior  to  any 
externally  applied  compressive  forces  or  the  beginning  of  a  subsequent 
load  cycle. 

Due  to  the  proposed  experimentation,  it  was  anticipated  that 
approximately  10-15  nodes  would  be  released  during  the  computer 
analysis,  thus  allowing  for  crack  growth  approximately  one  and  one-half 
to  two  crack  lengths .  The  release  of  s  node  was  assumed  to  occur  only 
after  a  steady  state  condition  was  reached  in  front  of  the  crack  tip. 
Since  it  was  not  feasible  to  attain  true  steady  state  conditions  prior 
to  the  release  of  each  node,  a  preliminary  study  was  accomplished  to 
determine  an  optimal  interval  between  subsequent  node  releases.  A 
single-edge  cracked  specimen  was  loaded  at  approximately  60  percent 
material  yield  strength  and  load  ratios  of  both  -1.0  and  0.1  were 
applied  in  two  separate  cases.  Crack  tip  stress  distributions  ahead  of 
the  crack  tip  and  displacements  behind  the  crack  tip  remained  virtually 
unchanged  after  one  complete  load  cycle  in  both  instances.  Plastic 
strains  ahead  of  the  crack  tip,  however,  required  several  cycles  to 
reach  a  near  steady  state  condition.  In  particular,  after  two  complete 
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load  cycles,  the  change  in  the  average  amount  of  plastic  strains  in  the 
three  elements  immediately  in  front  of  the  crack  tip  was  less  than  four 
percent  (Figure  3.4).  Additional  cycles  were  therefore  deemed 
unnecessary  when  considering  the  amount  of  crack  growth  anticipated. 

For  this  reason,  during  the  computer  analysis  of  the  single-edge  cracked 
specimens,  a  node  is  released  during  one  load  cycle  and  another  complete 
load  cycle  is  run  to  allow  for  the  stabilization  of  stress  and  strain 
conditions.  On  the  following  cycle,  data  are  collected  as  required  and 
the  procedure  (Figure  3.5)  is  repeated  until  the  desired  crack  growth  is 
attained . 

Finite  Element  Mesh 

In  consideration  of  the  proposed  experimentation  and  the  results  of 

preliminary  studies,  a  finite  element  mesh  (Figure  3.6)  consisting  of 

constant  strain  triangular  elements  was  developed.  Due  to  symmetry,  half 

of  the  specimen  has  been  modeled  and  necessary  displacement  boundary 

conditions  applied  at  the  nodes  along  the  axis  of  symmetry.  The  entire 

mesh  includes  389  elements  and  231  nodes.  The  fine  portion  of  the  mesh 

in  the  vicinity  of  the  crack  tip,  shown  in  Figure  3.7,  consists  of 
-4 

elements  1.6  x  10  inches  in  height  and  width  and  thus  an  area  of  1.28 

-  8  -  3 

x  10  inches  squared.  These  elements  extend  2.56  x  10  inches  ahead 

of  the  crack  tip  (approximately  two  and  one-half  crack  lengths)  in  order 

to  accomodate  the  experimental  crack  growth  data. 

The  size  of  the  fine  mesh  elements  is  based  upon  previous  work, 

experimental  measurement  accuracies,  one  -  dimensional  analyses,  and 

stress  intensity  factor  calculations.  Initial  estimates  were  obtained 

by  comparing  crack  tip  element  sizes  and  initial  crack  lengths  in 
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Figure  3.4 


Changes  in  Plastic  Strain 
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Figure  3.5 

Finite  Element  Crack  Growth  Procedure 
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Figure  3.6 

Finite  Element  Mesh  -  Single  Edge  Crack 
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Figure  3.7 

Crack  Tip  Finite  Element  Mesh 
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previous  studies.  Henkel  (68]  used  elements  measuring  7.8125  x  10 
inches  in  height  and  width  for  a  center  cracked  panel  (2a  -  .2734 
inches)  to  study  crack  closure  characteristics  under  cyclic  loading. 
Hinnerichs  [71]  employed  the  same  element  size  to  study  viscoplastic 
creep  in  several  center  cracked  panels  (2a  -  .2734,  .4733,  and  .6234 
inches) .  Comparative  elements  for  a  single-edge  cracked  specimen  (a  - 
.001  inch)  based  on  initial  crack  length  would  then  measure 
approximately  6.0  x  10  ^  inches.  In  addition,  data  obtained  from 
experimental  results  [22]  and  studies  of  the  methods  involved  indicated 
an  accuracy  of  lfxm  (3.9  x  10  inch)  in  measuring  crack  lengths  on  the 
order  of  25jnn  (.98  x  10  ^  inch)  [78]. 

A  one-dimensional  uniaxial  finite  element  analysis  was  also 
accomplished  using  small  elements  (2.0  x  10  ^  inch)  to  determine  if  the 
Visco  II  code  and  associated  computer  resources  could  accomodate  a 
small-sized  element.  In  question  was  the  ability  to  maintain  a 
sufficient  level  of  accuracy  in  the  solution  algorithms  when  element 
size  was  reduced.  Proposed  element  sizes  (in  terms  of  linear 
dimensions)  for  this  research  effort  were  over  1500  times  smaller  than 
any  element  previously  used.  By  applying  several  stress  values  under 
fully  reversed  cyclic  loads,  it  was  determined  that  the  size  of  an 
element  was  not  a  limiting  factor  as  long  as  the  resultant  stresses  did 
not  significantly  exceed  the  material  yield  strength.  This  indicates 
that  the  element  size  would  certainly  be  a  factor  when  the  final  mesh 
was  developed.  Since  the  stresses  in  a  constant  strain  triangle  are 
calculated  at  its  centroid,  as  the  element  size  is  reduced,  the  point  of 
stress  calculation  grows  closer  to  the  crack  tip.  This  results  in 
higher  crack  tip  element  stresses  for  a  given  external  load.  As  an 
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example,  in  an  elastic  analysis  the  stress  in  the  y-direction  in  the 
vicinity  of  the  crack  tip  is  related  by  [79,80]: 

oy  -  (KI/v/2jrr)cos0/2(l+sinfi/2sin30/2)  (3.5) 

For  a  single-edge  cracked  specimen,  Kj ,  the  stress  intensity  factor  can 
be  approximated  by  [79]: 

Kj  -  [1.99  -  .41(a/W)  +  18.7(a/W)2  -  38.48(a/W)3 

+  53.85(a/W)A]  aa1/2  (3.6) 

For  the  specimen  investigated  here,  this  reduces  to: 

K  -  . 062357609 cr  (3.7) 

Thus,  for  a  crack  tip  element  of  size  2.0  x  10  3  inch,  the  approximate 

stress  is  8. 472 O’  and  an  element  size  of  4.0  x  10  5  inch  has  an 

approximate  crack  tip  stress  of  5.99  O'  . 

As  a  result  of  the  above  considerations  and  the  requirement  for 

crack  growth  of  approximately  2-3  crack  lengths,  the  mesh  with  fine 

-4 

elements  measuring  1.6  x  10  inches  was  constructed.  Although  smaller 
element  sizes  may  ha'  e  been  utilized,  particularly  at  loads  of  60 
percent  yield  strength,  the  larger  load  of  90  percent  yield  strength 
precipitated  the  use  of  these  elements.  In  addition,  a  mesh  including 
smaller  elements  over  the  same  region  would  require  substantially  more 
elements  and  incur  excessive  computer  resources  to  simulate  crack  growch 
over  a  similar  distance. 

Finally,  an  external  load  was  applied  to  the  final  mesh  and  a 
calculation  of  the  stress  intensity  factor  performed.  For  an  externally 
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applied  stress  of  .04913  KSI,  R ,  was  analytically  determined  to  be  3.088 
1/2 

KSI-IN  using  equation  (3.7).  Using  an  extrapolation  technique  with 

1/2 

the  finite  element  mesh,  a  value  of  3.014  KSI-IN  was  found  (Figure 
3.8).  This  indicated  an  error  of  less  than  2.5  percent. 

Experimentation 

Specimens  composed  of  a  high  strength  titanium  alloy,  TI-6246,  and 
whose  geometry  is  shown  in  Figure  3.9,  were  subjected  to  cyclic  loading 
in  experiments  conducted  by  the  WRDC  Materials  Laboratory  [22,23]. 
Several  values  of  load  ratio,  R  (ratio  of  minimum  load  amplitude  to 
maximum  load  amplitude),  and  various  stress  levels  were  applied.  In 
particular,  load  ratios  of  -1.0  (fully  reversed  loading)  and  0.1  were 
investigated  in  conjunction  with  maximum  load  amplitudes  resulting  in 
nominal  stress  values  of  both  60  and  90  percent  material  yield  strength. 
Under  these  conditions,  growth  data  for  small  surface  flaws  were 
obtained.  The  particular  flaw  sizes  monitored  for  these  load  conditions 
are  shown  in  Table  3.1. 

The  analysis  of  the  experimental  data  involves  a  three-dimensional 
elastic  analysis  of  a  surface  flaw;  However,  it  is  proposed  that 
experimental  results  may  be  approximated  using  a  two-dimensional  finite 
element  analysis  of  a  single-edge  cracked  specimen.  The  crack  length, 
a,  in  this  case  would  be  equal  to  one-half  the  actual  surface  flaw 
length  (Figure  3.10).  Although  the  dimensionality  between  the  two  • 
cases  is  different,  the  finite  element  model  should  yield  comparable 
results  as  long  as  the  surface  flaw  is  semi-circular  (aspect  ratio  - 
a/r  -  1.00)  and  through- the  -  thickness  stresses  are  uniform. 

Experimental  results  confirmed  the  first  of  these  requirements  in  that 


MESH:  389  ELEMENTS 


Figure  3.8 

Stress  Intensity  Factor  Interpolation 
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Figure  3.9 


Figure  3.10 

Single -Edge  Cracked  Specimen 
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Figure  3.11 

Load-Time  Pattern  (R  -  -1.0) 


Figure  3.12 

Load-Time  Pattern  (R  -  0.1) 
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Table  3.1 


Experimental  Crack  Sizes 


LOAD 

R  - 

-1.0 

R 

-  0.1 

a.~  ,60a 

ys 

44  - 

530pm 

60  - 

1070pm 

.001732  - 

.02086  in 

.002362  - 

.04213 

in 

a  =  ,90a 
ys 

44  - 

860pm 

44  - 

814pm 

.001732  - 

.03386  in 

.001732  - 

.03205 

in 

Table  3.2 
Loading  Data 


-  ■ 

SPECIMEN 

a 

LOAD  (LB) 

KjrKSI-IN1/2) 

TI-6246 

.60  O’ 

10650.0 

6.69 

ys 

(106.5  KSI ) 

TI  -6246 

,90a 

15975.0 

10.04 

ys 

(159.75  KSI) 

Inconel  718 

.90a 

11700.0 

7.354 

ys 

(117.0  KSI) 
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the  surface  flaw  remained  nearly  semi-circular  (a/r  -  1.03)  and 

subsequent  analysis  by  the  Materials  Laboratory  verified  the  existence 

of  a  uniform  stress  distribution  through  the  specimen  thickness  [22]. 

Loads,  as  previously  stated,  Included  nominal  stress  values  of 

approximately  60  and  90  percent  yield  strength  in  conjunction  with  load 

ratios  of  -1.0  and  0.1.  For  TI-6246  (yield  strength  -  168  KSI),  the 

applied  loads  in  the  finite  element  analysis  are  shown  in  Table  3.2  and 

1/2 

resulted  initial  maximum  stress  intensity  factors  of  6.69  KSI-IN  and 

1/2 

10.04  KSI-IN  when  evaluated  using  equation  (3.6).  The  loads  were 
applied  in  a  triangular  saw-toothed  load- time  pattern  at  a  frequency  of 
1.0  Hz.  These  patterns,  for  the  two  load  ratios,  are  shown  in  Figures 
3.11  and  3.12.  For  Inconel  718  (yield  strength  -  130  KSI),  the  applied 
load  is  shown  in  Table  3.2.  A  nominal  stress  value  of  90  percent  yield 
strength  and  a  load  ratio  of  0.1  were  maintained,  whereas  the  load 
frequency  was  assumed  to  be  1.0  and  0.01  Hz. 
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IV.  Results  and  Discussion 


The  effects  of  plasticity  induced  closure  on  the  fatigue  growth  of 
short  cracks  was  investigated  by  analyzing  a  single-edge  cracked 
specimen  using  the  Visco  II  finite  element  code.  A  specimen  model 
consisting  of  TI-6246  at  room  temperature  was  subjected  to  four  separate 
load  conditions  in  which  the  load  ratio  and  nominal  stress  values  were 
varied  while  the  load  frequency  remained  constant  at  1.0  Hz.  A  second 
model  composed  of  high  temperature  Inconel  718  was  analyzed  using  two 
cases  where  the  load  ratio  and  nominal  stress  value  remained  constant, 
but  load  frequency  was  varied  in  order  to  observe  the  viscoplastic 
effects.  These  case  studies  are  outlined  in  Table  4.1.  In  each  of  these 
cases,  comparisons  are  made  between  the  case  in  which  a  crack  is  allowed 
to  grow  from  an  initial  crack  length  of  .001  inch  to  some  final  crack 
length  (.0026  in  for  TI-6246  and  .0018  in  for  Inconel  718),  and  the  case 


Table  4.1 

Finite  Element  Case  Studies 


Case 

Material 

Stress 

R 

af 

Freq 

<a/V 

(in) 

(Hz) 

A 

TI-6246 

.60 

-1.0 

.0026 

1.0 

B 

TI-6246 

.60 

0.1 

.0026 

1.0 

C 

TI-6246 

.90 

-1.0 

.0026 

1.0 

D 

TI-6246 

.90 

0.1 

.0026 

1.0 

E 

Inconel  718 

.90 

0.1 

.0018 

1.0 

F 

Inconel  718 

.90 

0.1 

.0018 

.01 
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where  the  initial  crack  length  is  set  equal  to  that  of  the  grown  crack 
in  an  unde  formed  material  where  crack  extension  has  not  occurred.  A 
graphic  representation  of  this  finite  element  crack  growth  is  shown  in 
Figure  4.1.  Crack  growth  is  obtained  through  the  use  of  the  node 
release  algorithm  discussed  in  Chapter  2.  By  analyzing  the  two  cracks 
in  this  manner,  the  effects  of  fatigue  crack  growth  history  and  the 
formation  of  a  plastic  wake  can  be  better  isolated. 

For  the  two  materials,  data  have  been  obtained  from  the  finite 
element  analyses  and  are  presented  in  several  categories.  These  include 
displacements  along  the  crack  surface,  stress-strain  values  immediately 
in  front  of  the  crack  tip,  stress  distributions  ahead  of  the  crack  tip, 
plastic  strains  along,  as  well  as  in  front  of,  the  crack,  closure  data, 
and  load-displacement  relationships. 

Results  for  Tl-6246 

The  single-edge  cracked  specimen  model  consisting  of  TI-6246,  a 

high  strength  titanium  alloy  at  room  temperature,  was  numerically 

subjected  to  external  loads  of  10650  and  15975  pounds  resulting  in 

nominal  stress  values  of  approximately  60  and  90  percent  of  the  material 

yield  strength,  respectively.  At  the  initial  crack  length  of  .001  inch, 

these  loads  represent  stress  intensity  factors  of  approximately  6.69  and 

1/2 

10.04  KSI-IN  In  all  but  one  case,  ten  finite  element  nodes  were 

released  during  the  analysis  which  resulted  in  a  crack  extension  of  • 
.0016  inch  (a^.  -  .0026  inch).  In  case  B,  a  total  of  fifteen  nodes  were 
released  (a^  =  .00324  inch).  Results  and  trends  observed  at  this  last 
crack  length,  however,  did  not  differ  significantly  from  those  at  the 
shorter  final  crack  length.  Therefore,  the  results  obtained  from  all 
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four  cases  are  presented  for  the  same  final  crack  length  (a^.  -  .0026 
inch)  for  reasons  of  comparison  and  continuity. 

The  y-displacement  profiles  behind  the  crack  tip  at  minimum 
loading  show  significant  differences  between  the  propagated  and 
non-propagated  cracks.  In  cases  A  and  C,  Figures  4.2  and  4.4,  where 
the  load  ratio  is  -1.0,  the  displacements  in  the  non-propagated  crack 
are  much  greater  and  occur  immediately  behind  the  crack  tip.  The 
displacements  in  the  propagated  cracks  occur  at  a  distance  of  .00176  in 
behind  the  crack  tip  (in  the  vicinity  of  the  original  crack  tip)  and 
are  considerably  less.  Similar  results  are  shown  in  Figures  4.3  and 
4.5  for  the  cases  where  the  load  ratio  is  0.1.  In  these  cases,  the 
non-propagated  crack  remains  fully  open  at  minimum  load  and  the  maximum 
displacement  in  the  vicinity  of  the  crack  tip  occurs  at  the  node 
directly  behind  the  tip.  The  displacement  profile  of  the  propagated 
crack  indicates  that  the  crack  remains  closed  for  several  nodes 
adjacent  to  the  crack  tip  and  then  opens  gradually  up  to  the  location  of 
the  original  tip.  At  this  point,  the  displacement  increases 
significantly.  These  findings  show  that  for  the  non-propagated  crack, 
the  crack  tip  is  fully  open  and  maximum  displacement  occurs  at  the  node 
directly  behind  the  crack  tip.  For  the  propagated  crack,  however,  the 
crack  remains  closed  for  several  nodes  behind  the  crack  tip  and  then 
opens  gradually  as  distance  from  the  crack  tip  increases.  The 
discontinuity  which  occurs  in  the  vicinity  of  the  original  crack  tip 
(Figures  4.3  and  4.5)  indicates  that  under  cyclic  loading  the 
displacements  take  on  unusually  high  values  directly  behind  the  crack 
tip  when  no  history  from  prior  fatigue  cycling  and  crack  extension  are 
considered.  It  is  unrealistic  to  represent  displacement  profiles  behind 
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Af  =  .00260  in 
MINIMUM  LOAD 


Fig  4.2:  Displacements  Behind  Crack  Tip  (Case  A/Min  Load) 
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Af  =  .00260  in 
MINIMUM  LOAD 


Fig  4.3.  Displacements  Behind  Crack  Tip  (Case  B/Min  Load) 
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Fig  4  4:  Displacements  Behind  Crack  Tip  (Case  C/Min  Load) 
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DISTANCE  (10“3  IN) 


Fig  4.5:  Displacements  Behind  Crack  Tip  (Case  D/Min  Load) 
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a  propagating  crack  in  a  finite  element  model  which  does  not  include  the 
cumulative  effects  of  crack  propagation. 

Analysis  of  the  y-displacement  profiles  behind  the  crack  tip  at 
maximum  loading  yields  similar  results  to  those  observed  at  minimum 
loading.  In  each  case,  Figures  4.6  to  4.9,  the  displacements  directly 
behind  the  crack  tip  observed  in  the  non-propagated  cracks  are  greater 
than  those  of  the  propagated  cracks.  The  existence  of  a  plastic  wake 
resulting  from  cyclic  loading  and  crack  extension  reduces  the  amount  of 
displacement  along  the  crack  surface.  The  displacements  which  are 
measured  past  the  original  crack  tip,  i.e.,  further  than  .0016  in,  are 
nearly  identical  for  all  four  cases.  The  effects  of  cyclic  loading  and 
crack  extension  do  not  occur  over  the  original  .001  inch  crack  length. 
The  larger  displacements  always  correspond  to  the  non-propagated  crack. 
The  difference  between  these  profiles  corresponds  to  the  plastic  wake 
which  forms  behind  the  propagating  crack.  Therefore,  the  plastic  wake 
is  not  modeled  unless  nodes  are  released  in  the  finite  element  model 
simulating  crack  growth. 

The  stress/strain  curves  for  elements  immediately  in  front  of  the 
crack  tip  are  displayed  in  Figures  4.10  through  4.13.  In  each  case, 
the  maximum  and  minimum  stresses  observed  do  not  differ  significantly 
between  the  propagated  and  non-propagated  cracks.  The  stresses, 
therefore,  do  not  appear  to  be  affected  by  crack  growth  history.  The 
strains  exhibited  by  the  non-propagated  cracks,  however,  are  greater  in 
each  case.  Under  a  nominal  stress  value  of  60  percent  yield  strength, 
Figures  4.10  and  4.11,  the  non-propagated  crack  exhibits  strains 
approximately  25  percent  greater  than  the  propagated  crack.  At  higher 
loads,  as  in  Figure  4.13,  this  increase  in  strain  is  nearly  doubled. 
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Fig  4.6:  Displacements  Behind  Crack  Tip  (Case  A/Max  Load) 
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Fig  4.7:  Displacements  Behind  Crack  Tip  (Case  B/Max  Load) 
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Fig  4.8:  Displacements  Behind  Crack  Tip  (Case  C/Max  Load) 
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Fig  4.9:  Displacements  Behind  Crack  Tip  (Case  D/Max  Load) 
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Fig  4.10:  Stress/Strain  Ahead  of  Crack  Tip  (Case  A) 
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Fig  4.11:  Stress/Strain  Ahead  of  Crack  Tip  (Case  B) 
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The  strains  immediately  ahead  of  the  crack  tip  are  therefore  related  to 
crack  growth  history.  The  plastic  wake  formed  as  a  result  of  crack 
growth  limits  the  amount  of  strain  exhibited.  Since  the  non-propagated 
crack  is  not  subjected  to  this  limitation,  it  shows  higher  strain 
values . 

As  illustrated  in  Figures  4.10  through  4.13,  the  stresses  ahead  of 
the  crack  tip  do  not  differ  significantly  between  propagated  and 
non-propagated  cracks.  Figures  4.14  through  4.17  show  the  stress 
profiles  ahead  of  the  crack  tip.  The  y-stress  at  full  positive  load  is 
greatest  near  the  crack  tip  and  reduces  gradually  as  the  distance  from 
the  crack  tip  is  increased.  Ac  minimum  loading,  the  negative  stresses 
increase  rapidly  over  the  first  five  elements  in  front  of  the  crack  tip 
and  then  level  off  for  each  case  except  for  case  C  where  the  load  ratio 
is  -1.0  and  nominal  stress  is  .90  yield  strength.  Here  the  high 
compressive  forces  are  nearly  constant  over  the  first  six  elements  and 
then  begin  to  increase  gradually. 

The  distribution  of  plastic  strain  in  the  y-direction  at  maximum 
tensile  load  for  the  four  cases  is  presented  in  F’fures  4.18  to  4.21. 

The  differences  between  the  propagated  and  non-propagated  cracks  are 
obvious.  In  each  instance,  the  non-propagated  crack  exhibits  no  plastic 
strain  along  the  crack  length  and  has  a  slightly  greater  strain  at  the 
crack  tip.  The  non-propagated  cases  show  an  increase  in  plastic  strain 
as  the  crack  tip  is  approached  and  do  not  exhibit  as  great  a  strain  in 
the  crack  tip  vicinity.  As  the  crack  is  allowed  to  grow,  plastic 
strains  are  accumulated  along  the  crack  surface.  These  residual  strains 
restrict  the  amount  of  plastic  strain  observed  in  the  crack  tip  element. 
Ahead  of  the  crack  tip,  the  elastic  strains  in  tne  non-prepagated  cracks 
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Fig  4.14:  Stresses  Ahead  of  Crack  Tip  (Case  A) 
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Fig  4.15:  Stresses  Ahead  of  Crack  Tip  (Case  B) 
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Fig  4.17:  Stresses  Ahead  of  Crack  Tip  (Case  D) 


.60  YIELD  STRENGTH 


CD 


Fig  4.18:  Plastic  Strain  Along  X-Axis  (Case  A) 
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Fig  4.19:  Plastic  Strain  Along  X-Axis  (Case  B) 
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Fig  4.21:  Plastic  Strain  Along  X-Axis  (Case  D) 
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decrease  asymptotically  and  approach  those  of  the  propagated  cracks 
within  several  elements.  In  addition,  the  amount  of  plastic  strain 
exhibited  appears  to  be  a  function  of  both  the  nominal  stress  applied 
and  the  load  ratio.  In  ca$es  B  and  D,  in  Figures  4.19  and  4.21  where 
the  load  ratio  is  -1.0,  the  nominal  stress  is  increased  from  60  to  90 
percent  yield  strength.  The  plastic  strain  observed  in  the  propagated 
cracks  rises  dramatically  in  both  cases.  The  plastic  strain  then 
appears  to  vary  directly  with  load  magnitude.  When  the  load  ratio  is 
increased,  a  similar  condition  exists.  At  a  load  ratio  of  0.1,  the 
plastic  strains  along  the  crack  length  are  greater  than  those  with  a 
load  ratio  of  -1.0.  The  compressive  forces  present  with  the  negative 
load  ratio  limit  the  accumulation  of  plastic  strain. 

Analysis  of  the  y-direction  plastic  strains  at  maximum  load  yield 
several  findings.  No  residual  strains  are  exhibited  behind  the  crack 
tip  for  the  non-propagated  crack.  Ahead  of  the  crack  tip,  however, 
these  strains  slightly  exceed  those  of  the  propagated  crack  near  the  tip 
and  then  asymptotically  approach  those  of  the  propagated  crack  at 
greater  distances  from  the  crack  tip.  The  distance  over  which  the 
plastic  strains  are  measureable  ahead  of  the  crack  tip,  i.e.,  the  size 
of  the  plastic  zone,  is  approximately  equal  to  the  initial  crack  length 
of  .001  inch.  This  condition  is  typical  of  short  crack  phenomena.  The 
plastic  strains  immediately  ahead  of  the  crack  tip  are  higher  for  the 
non-propagated  crack  than  for  the  propagated  crack.  Further,  by 
observing  the  strains  behind  the  crack  tip,  it  is  seen  that  the  plastic 
residual  strains  increase  significantly  as  the  nominal  applied  stress 
increases  from  0.6  to  0.9  of  the  yield  stress  in  the  case  of  R  -  0.1. 
However,  for  fully  reversed  loading,  R  -  -1.0,  the  residual  plastic 
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strains  are  lower  and  become  smallest  at  the  highest  stress  level.  This 
indicates  that  severe  compression  develops  in  the  wake  under  fully 
reversed  loading  and  the  plastic  wake  is  compressed  the  most  at  the 
higher  stress  level  in  the  region  directly  behind  the  crack  tip. 

by  comparing  the  y-displacements  behind  the  crack  tip  with 
corresponding  residual  strain  data,  several  findings  are  made.  The 
differences  in  the  y-displacements  for  the  non-propagated  crack 
(larger  displacements)  and  the  propagated  crack  (smaller  displacements) 
correspond  to  the  plastic  wakes  behind  the  crack  tip  formed  during  crack 
growth.  For  the  lower  stress  level,  i.e.,  0.6  yield  stress,  the  plastic 
wake  is  larger  for  R  -  0.1  than  for  the  case  where  R  -  -1.0.  This 
agrees  with  the  amount  of  residual  plastic  strain  observed  in  those 
cases.  At  the  higher  stress  level,  i.e.,  0.9  yield  stress,  the  larger 
plastic  wake  at  R  -  0.1  is  observed  in  both  the  displacement  profile  and 
the  strain  profile.  For  R  -  -1.0,  larger  displacements  occur  at  maximum 
load  but  the  size  of  the  plastic  wake  is  similar  to  those  of  -  ,e  lower 
stressvalues  near  the  crack  tip.  At  larger  distances  from  th<-  crack 
tip,  residual  strains  decrease  while  maximum  displacements  become  larger 
for  the  propagated  crack  than  for  the  non-propagated  crack.  The  highest 
amount  of  plastic  deformation  occurs  for  the  higher  stress  level  (0.9 
yield  stress)  at  R  -  0.1. 

Evidence  of  the  increased  amount  of  plasticity  along  the  crack 
length  due  to  fatigue  crack  growth  is  observed  in  several  forms  of 
closure  data.  Figures  .22  through  A. 25  present  crack  closure  and 
opening  data  as  a  function  of  crack  length.  Here,  P(clo)  is  the  load 
at  which  some  point  on  the  crack  closes.  In  the  finite  element 
analysis,  this  occurs  when  a  node  on  the  crack  surface  exhibits  a 
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Fig  4.22:  Crack  Closure  (Case  A) 
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Fig  4.23:  Crack  Closure  (Case  B) 
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Fig  4.24:  Crack  Closure  (Case  C) 
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Fig  4.25:  Crack  Closure  (Case  D) 
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negative  displacement,  thereby  initiating  the  crack  closure  algorithm. 
P(ope)  is  the  first  point  during  reloading  when  all  nodes  along 
the  crack  surface  exhibit  positive  displacements.  At  60  percent 
material  yield  strength,  both  closure  and  opening  values  continue  to 
increase  as  the  crack  propagates.  As  the  crack  grows  and  the  plastic 
wake  forms,  there  is  a  corresponding  increase  in  residual  compressive 
forces  within  the  material  during  unloading.  As  a  result,  less  external 
compression  is  required  to  close  the  crack  and  greater  tension  is 
required  before  the  crack  reopens  on  subsequent  loading.  In  cases  C  and 
D,  Figures  4.24  and  4.25,  the  nominal  stresses  are  greater.  The  closure 
data  here  does  not  show  the  increases  that  are  evident  at  lower  stress 
values.  In  these  cases,  the  closure  and  opening  values  are  nearly 
constant.  There  appears  to  be  a  limit  in  closure  loads  which  is  a 
function  of  the  load  conditions.  For  fully  reversed  loading  as  in  case 
C,  the  closure  load  decreases  as  the  crack  extends.  As  the  size  of  the 
plastic  area  is  reduced  because  of  high  compressive  stresses,  the 
residual  compressive  stresses  are  also  reduced.  As  a  result,  greater 
forces  are  required  to  close  the  crack. 

The  difference  between  crack  opening  and  crack  closure  discussed 
above  is  most  noticeable  for  the  cases  where  R  -  -1.0  which  exhibit 
large  hysteresis  loops  under  fully  reversed  loading.  This  hysteresis  is 
largest  at  the  higher  stress  level  of  0.9  yield  stress  (Figures  4.10  and 
4.12).  Further,  closure  continues  to  develop  over  the  total  amount  of 
crack  extension  at  the  lower  stress  level  of  0.6  yield  stress.  If  the 
crack  does  not  propagate,  the  closure  load  at  the  final  crack  length  is 
approximately  the  same  as  for  the  crack  before  it  was  allowed  to  grow. 

At  higher  stress  levels  the  closure  load  does  not  vary  considerably  as 
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the  crack  extends.  At  the  higher  stress  levels  the  plastic  deformation 
occurring  behind  the  crack  tip  reaches  an  equilibrium  condition  with 
little  or  no  extension.  This  is  also  observed  for  the  case  where  the 
nominal  stress  is  0.9  yield  stress  and  fully  reversed  loading  occurs. 

The  closure  load  for  the  longer  crack  length  with  no  plastic  strain 
history  is  nearly  the  same  for  the  propagated  crack.  At  this  higher 
stress  level,  but  where  R  -  0.1,  the  closure  load  reaches  a  value 
slightly  greater  than  the  minimum  load  after  crack  extension  and  then 
remains  constant.  For  the  non-propagated  crack,  no  closure  develops  and 
therefore  the  closure  load  must  be  less  than  the  minimum  load  applied. 

Another  method  of  measuring  closure  is  to  observe  the 
closure  load  as  a  function  of  crack  length.  Closure  load  has  been 
converted  to  closure  stress  intensity  using  the  elastic  stress  intensity 
solution  for  the  single  edge  cracked  specimen  and  these  data  have  been 
plotted  against  K(max)  for  all  four  cases  in  Figures  4.26  to  4.29.  At 
lower  applied  stresses  (Figures  4.26  and  4.27),  K(clo)  and  K(ope)  for 
the  propagating  crack  increase  with  K(max) .  At  higher  applied  stresses 
(Figures  4.28  and  4.29),  the  same  is  not  true.  At  a  load  ratio  of  -1.0, 
the  values  of  K(clo)  and  K(ope)  decrease  with  an  increase  in  K(max)  and 
remain  nearly  constant  at  a  load  ratio  of  0.1.  The  higher  compressive 
stresses  ahead  of  the  crack  tip  result  in  lower  values  of  K(clo)  and 
K(ope)  . 

For  plasticity  induced  closure,  analytical  and  experimental  results 
indicate  that  K(clo)  or  K(ope)  increases  linearly  as  K(max)  increases 
[32],  For  R  -  0.1,  K(clo)  increases  with  K(max)  for  the  lower  applied 
stress.  For  the  higher  applied  stress,  K(clo)  decreases  with  increases 
in  K(max) .  This  indicates  that  there  is  no  single-valued  function  of 
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Fig  4.26:  Crack  Closure  As  a  FVinction  of  K  (Case  A) 
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Fig  4.27:  Crack  Closure  As  a  Function  of  K  (Case  B) 
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Fig  4.28.  Crack  Closure  As  a  Function  of  K  (Case  C) 
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Fig  4  29:  Crack  Closure  As  a  Function  of  K  (Case  D) 
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K(clo)  against  K(max)  for  R  -  0.1.  At  the  higher  stress  levels  there 
are  compressive  stresses  generated  ahead  of  the  crack  tip  because  of  the 
higher  stress  range  coupled  with  the  maximum  stress  being  limited  by  the 
yield  stress.  This  results  in  a  residual  deformation  pattern  caused  by 
cyclic  plasticity  leading  to  lower  values  of  K(clo) .  At  a  stress  ratio 
of  R  -  -1.0  and  at  the  higher  stress  level,  the  value  of  K(ope)  remains 
relatively  constant  at  a  sightly  negative  value  while  K(clo)  decreases 
with  increasing  crack  extension.  K(ope)  and  K(clo)  for  the  non¬ 
propagating  crack  do  not  differ  significantly  from  the  propagating  crack 
at  the  final  crack  length.  Closure  here  appears  to  be  a  result  of 
cyclic  plastic  deformation  at  or  near  the  crack  tip  rather  than  as  a 
result  of  the  residual  plastic  wake  behind  the  propagating  crack.  At 
the  lower  stress  level,  however,  K(ope)  and  K(clo)  for  the  non¬ 
propagating  crack  differ  significantly  from  the  propagating  crack.  Here 
K(ope)  increases  faster  than  K(clo)  as  K(max)  increases. 

As  previously  discussed,  closure  may  be  analytically  determined  as 
the  point  where  nodal  displacements  become  either  negative  or  zero  in 
the  finite  element  code.  Experimentally,  however,  closure  data  is 
normally  determined  through  the  use  of  load-displacement  curves  obtained 
from  measurements  taken  from  a  test  specimen.  Out  of  necessity,  these 
measurements  are  usually  taken  at  or  near  the  crack  mouth.  Figures  4.30 
through  4.53  present  load-displacement  and  load-delta  displacement  for 
three  nodes  along  the  crack  surface.  Node  192  represents  the  crack 
mouth,  whereas  nodes  205  and  206  are  those  immediately  adjacent  to  the 
crack  tip.  Figure  3.7  shows  this  node  numbering  scheme.  The  closure 
load  determined  from  load- displacement  data  is  that  point  at  which  the 
curve  becomes  non-linear.  At  points  distant  from  the  crack  tip,  this 
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Fig  4.30:  Load  -  Displacement  (Case  A/Elem  192) 
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Fig  4.31:  Load  -  Delta  Displacement  (Case  A/Elem  192) 


88 


Fig  4.32:  Load  -  Displacement  (Case  A/Elem  205) 
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Fig  4.33:  Load  -  Delta  Displacement  (Case  A/Elem  205) 
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Fig  4.34:  Load  -  Displacement  (Case  A/Elem  206) 
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Fig  4.35:  Load  -  Delta  Displacement  (Case  A/Elem  206) 
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Fig  4.36:  Load  -  Displacement  (Case  B/Elem  192) 
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Fig  4.37:  Load  -  Delta  Displacement  (Case  B/Elem  192) 
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Fig  4  38:  Load  —  Displacement  (Case  B/Elem  205) 
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Fig  4.39:  Load  -  Delta  Displacement  (Case  B/Elem  205) 
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Fig  4.40:  Load  -  Displacement  (Case  B/Elem  206) 
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Fig  4.41:  Load  -  Delta  Displacement  (Case  B/Elem  206) 
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Fig  4.42:  Load  -  Displacement  (Case  C/Elem  192) 
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Fig  4.43:  Load  -  Delta  Displacement  (Case  C/Elem  192) 
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Fig  4.44:  Load  -  Displacement  (Case  C/Elem  205) 
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Fig  4.45:  Load  -  Delta  Displacement  (Case  C/Elem  205) 
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Fig  4.46:  Load  -  Displacement  (Case  C/Elem  206) 
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Fig  4.47:  Load  -  Delta  Displacement  (Case  C/Elem  206) 
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Fig  4.48.  Load  -  Displacement  (Case  D/Elem  192) 
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Fig  4.49:  Load  -  Delta  Displacement  (Case  D/Elem  192) 
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Fig  4.50:  Load  -  Displacement  (Case  D/Elem  205) 
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Fig  4.51:  Load  -  Delta  Displacement  (Case  D/Elem  205) 
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Fig  4.52:  Load  -  Displacement  (Case  D/Elem  206) 
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point  is  difficult  to  determine  because  of  the  accumulation  of 
plasticity  effects  between  the  crack  tip  and  the  point  of  displacement 
measurement.  In  each  of  the  four  cases  presented,  it  is  observed  that 
the  nearer  to  the  crack  tip  that  measurements  are  taken,  the  more 
distinct  the  change  in  linearity  of  the  load-displacement  curve.  At 
nodes  immediately  behind  the  crack  tip,  very  accurate  determinations  of 
crack  closure  can  be  obtained.  Those  further  from  the  crack  tip  are 
less  reliable. 

An  improvement  to  the  load-displacement  method  of  experimentally 
determining  closure  is  to  modify  the  data  to  obtain  a  load-delta 
displacement  curve.  As  a  specimen  is  initially  unloaded,  a  measurement 
of  the  displacement  curve  is  obtained.  Deviations  from  the  analytically 
derived  slope  are  then  used  to  determine  when  closure  occurs.  Although 
this  technique  improves  the  accuracy  in  determining  the  point  of 
closure,  far- field  measurements  are  still  obscured  by  the  accumulation 
of  plasticity.  As  measurements  are  taken  closure  to  the  crack  tip,  more 
accurate  determinations  of  crack  closure  and  opening  are  obtained. 

Other  observations  regarding  crack  opening  and  closure  can  be  made 
by  further  analysis  of  the  load-delta  displacement  curves.  During  the 
time  the  crack  is  closed  and  external  load  is  increasing,  the  slope  of 
the  curve  remains  constant  until  the  crack  opens.  At  this  point,  there 
is  a  distinct  change  in  slope  and  the  value  of  P(ope)  is  easily 
obtained.  As  the  load  continues  to  increase,  the  slope  of  the 
load-delta  displacement  curve  is  again  constant  representing  the 
occurrence  of  elastic  strains.  When  the  applied  load  reaches  its 
maximum  value,  plastic  strains  begin  to  occur  and  a  third  change  in 
slope  is  observed. 
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During  the  unloading  portion  of  the  load  cycle,  the  specimen 
initially  undergoes  an  elastic  reduction  in  strain  which  is  represented 
by  the  vertical  slope  of  the  load-delta  displacement  curve.  At  the 
load  corresponding  to  P(ope),  a  portion  of  that  plastic  strain  is 
recovered  until  crack  closure  occurs.  This  is  indicated  by  a  reduction 
in  the  curve's  slope.  At  the  point  of  closure,  the  slope  is  once  more 
distinctly  reduced.  Furthermore,  the  difference  between  the  values  of 
P(clo)  and  P(ope) ,  coupled  with  the  amount  of  recoverable  plastic 
strain,  determine  the  extent  of  hysteresis  exhibited  in  the  load-delta 
displacement  curves. 

Results  for  Inconel  718 

The  single-edge  cracked  specimen  composed  of  high  temperature 
Inconel  718  was  subjected  to  a  constant  nominal  stress  value  of  90 
percent  of  the  material  yield  strength  and  a  load  ratio  of  0.1,  whereas 
the  frequency  was  varied  between  1.0  and  .01  Hz.  In  these  cases,  one 
crack  was  allowed  to  grow  from  an  initial  crack  length  of  .001  inch  to  a 
final  length  of  .0018  inch  (five  finite  element  nodes  released),  and 
the  second  crack  had  an  initial  crack  length  of  .0018  inch.  Through 
the  use  of  the  two  load  frequencies,  the  effects  of  viscoplasticity  on 
plasticity  induced  closure,  the  formation  of  a  plastic  wake,  and 
fatigue  crack  growth  were  investigated. 

Figures  4.54  through  4.57  show  the  displacements  behind  the  crack 
tip  at  both  minimum  and  maximum  loads.  As  expected,  the  displacements 
were  less  for  the  propagated  cracks  over  the  crack  growth  region.  The 
differences  in  the  displacements  due  to  the  nature  of  propagating  and 
non-propagating  cracks  are  nearly  the  same  for  each  frequency.  The 
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Fig  4.54:  Displacements  Behind  Crack  Tip  (Case  E/Min  Load) 
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Fig  4.55.  Displacements  Behind  Crack  Tip  (Case  I*/Min  Load) 
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Fig  4.56:  Displacements  Behind  Crack  Tip  (Case  E/Max  Load) 
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Fig  4.57:  Displacements  Behind  Crack  Tip  (Case  F/Max  Load) 
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magnitudes  of  the  displacements  observed  at  the  lower  frequency  (case 
F) ,  however,  are  greater  than  the  higher  frequency.  The  viscoplastic 
properties  at  various  load  frequencies  do  not  significantly  affect  the 
differences  in  propagated  and  non-propagated  crack  for  this  material. 

The  stress/strain  hysteresis  loops  shown  in  Figures  4.58  and  4.59 
show  that  stresses  ahead  of  the  crack  tip  are  not  significantly 
affected  by  previous  crack  growth.  Further,  these  stresses  do  not 
differ  with  changes  in  load  frequency.  The  strain  values  observed, 
however,  are  reduced  as  a  result  of  crack  propagation  and  the  formation 
of  the  plastic  wake.  The  amount  of  reduction  in  strain  resulting  from 
crack  growth  is  less  in  case  E,  indicating  that  the  higher  frequency 
does  not  allow  as  much  straining  as  the  lower  frequency. 

As  observed  in  the  hysteresis  loops  previously  discussed,  the 
stress  values  ahead  of  the  crack  tip  are  only  slightly  affected  by 
crack  growth  and  viscoplasticity.  In  Figures  4.60  and  4.61  the 
distribution  of  stresses  at  minimum  and  maximum  load  values  show  a 
slight  difference  due  to  the  fatigue  crack  growth  and  formation  of  the 
plastic  wake.  There  is  also  a  slight  increase  in  all  stress  values  in 
the  vicinity  of  the  crack  tip  for  the  higher  load  frequency  which 
appears  to  be  the  result  of  viscoplastic  material  behavior  of  IN-718  at 
high  temperature. 

Figures  4.62  and  4.63  display  the  accumulation  of  plastic  strain 
along  and  in  front  of  the  crack  tip.  As  the  crack  propagates,  the 
accumulation  of  plastic  strain  along  the  region  of  crack  growth  is 
evident  as  is  the  reduced  strain  immediately  in  front  of  the  crack  tip 
for  the  propagated  crack.  This  phenomenon  was  also  exhibited  in  the 
stress/strain  curves  previously  discussed  and  appears  to  be  a  result  of 
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Fig  4.58:  Stress/Strain  Ahead  of  Crack  Tip  (Case  E) 
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Fig  4.59.  Stress/Strain  Ahead  of  Crack  Up  (Case  F) 
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Fig  4.60:  Stresses  Ahead  of  Crack  Tip  (Case  E) 
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Fig  4.61:  Stresses  Ahead  of  Crack  Tip  (Case  F) 
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Fig  4.63:  Plastic  Strain  Along  X-Axis  (Case  F) 
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the  plastic  wake  behind  the  crack  tip.  Further,  by  superimposing  the 
plastic  strains  for  these  two  cases  as  in  Figure  4.64  the  ef'ect  of 
viscoplasticity  can  be  observed.  At  higher  load  frequencies  the 
strains  are  slightly  reduced  along  the  entire  crack  length  for  both 
propagated  and  non-propagated  cracks. 

Closure  data  as  related  to  crack  length  and  maximum  stress 
intensity  factor  are  shown  in  Figures  4.65  to  4.68.  In  each  case  the 
closure  and  opening  loads  are  reduced  as  the  crack  extends  and  K(max) 
increases.  The  viscoplastic  behavior  of  the  Inconel  718  material 
permits  the  growth  of  plastically  deformed  material  which,  in  turn, 
increases  residual  compressive  stresses.  As  a  result,  the  cracks  close 
at  a  higher  external  load  during  unloading  than  they  would  have  in  the 
absence  of  any  plastic  deformation.  Similarly,  during  the  subsequent 
loading  portion  of  the  fatigue  cycle,  crack  opening  occurs  at  a  higher 
external  load. 

Load-displacement  and  load-delta  displacement  curves  are  presented 
for  the  Inconel  718  specimens  in  Figures  4.69  through  4.80.  These  have 
been  obtained  for  three  node  locations.  The  crack  mouth  is  represented 
by  node  192,  and  the  area  immediately  behind  the  crack  tip  by  nodes  200 
and  201.  As  discussed  previously,  it  is  apparent  that  the  point  of 
closure  or  opening  is  barely  distinguishable  when  displacements  are 
taken  far-field,  such  as  at  the  crack  mouth.  As  the  location  of 
measurements  moves  closer  to  the  crack  tip  more  accurate  estimations  of 
crack  closure  and  opening  can  be  obtained. 

In  these  curves  for  the  Inconel  718  specimens,  the  same 
characteristics  as  were  evident  in  the  TI-6246  specimen  are  observed. 
During  loading,  the  load-delta  displacement  slope  remains  constant  as 
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Fig  4.64:  Affect  of  Frequency  on  Plastic  Strain 
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Fig  4.65:  Crack  Closure  (Case  E) 
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Fig  4.66:  Crack  Closure  (Case  F) 
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Fig  4.67:  Crack  Closure  As  a  Function  of  K  (Case  E) 
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Fig  4.68:  Crack  Closure  As  a  Function  of  K  (Case  F) 
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Fig  4.69:  Load  -  Displacement  (Case  E/Elem  192) 
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Fig  4.70:  Load  -  Delia  Displacement  (Case  E/Elem  192) 


131 


[  I  I  I  I  I 


oo  <o  c\}  o 

©  ©  ©  ©  - 

XVWd/d 

Fig  471:  Load  -  Displacement  (Case  E/Elem  200) 
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Fig  472:  Load  -  Delta  Displacement  (Case  E/Elem  200) 
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Fig  4.73:  Load  -  Displacement  (Case  E/Elem  201) 


134 


P/PMAX 


Fig  4.74:  Load  -  Delia  Displacement  (Case  E/Elem  201) 
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Fig  475:  Load  -  Displacement  (Case  F/Elem  192) 
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Fig  4.76:  Load  -  Delia  Displacement  (Case  F/Elem  192) 
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Fig  4.77:  Load  -  Displacement  (Case  I^/Elem  200) 
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Fig  4.78:  Load  -  Delta  Displacement  (Case  F/Elem  200) 
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Fig  4  79  Load  -  Displacement  (Case  F/Elem  201) 
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Fig  4.80:  Load  -  Delia  Displacement  (Case  £/Elem  201) 
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the  crack  remains  closed,  but  at  the  point  of  crack  opening,  it 
increases  dramatically.  From  this  point,  the  elastic  strains 
predominantly  occur  until  plastic  strains  become  evident  near  the  point 
of  maximum  load.  During  subsequent  unloading,  a  reduction  in  Lite 
elastic  strains  initially  occurs,  but  then  the  recovery  of  plastic 
strains  predominates  as  P(clo)  is  approached.  The  point  of  crack 
closure  is  indicated  by  a  sharp  change  in  the  slope  of  the  load-delta 
displacement  curve. 

The  effects  of  viscoplastic  material  properties  in  Inconel  718  at 
high  temperatures  can  be  further  seen  in  these  curves.  As  an  example, 
the  load-delta  displacement  curve  in  Figure  4.74  shows  that  the  width  of 
the  hysteresis  loop  is  approximately  .6  x  10  ^  inches.  This  represents 
measurements  taken  at  the  node  immediately  behind  the  crack  tip  with  a 
load  frequency  of  1.0  Hz.  Similar  measurements  at  a  cyclic  frequency  of 
.01  Hz  in  Figure  4.80  reflect  a  width  of  .85  x  10  6  inches.  A  greater 
amount  of  recoverable  plastic  strain  occurs  at  the  lower  frequency. 

Discussion 

Regardless  of  the  materials  involved,  here  TI-6246  at  room 
temperature  and  Inconel  718  at  1200°F,  as  a  crack  propagates  under 
cyclic  loading,  a  plastic  wake  is  formed  along  the  crack  surface.  This 
wake  is  composed  of  plastically  deformed  material  and  has  a  definite 
effect  of  material  behavior  in  the  vicinity  of  the  crack  tip.  In  order 
to  accurately  describe  this  behavior,  finite  element  analyses  must  be 
able  to  incorporate  the  development  and  effects  of  a  plastic  wake.  As 
seen  in  this  research,  finite  element  modeling  of  plasticity  induced 
closure  in  short  cracks  under  fatigue  loading  requires  some  amount  of 
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crack  extension  to  develop  closure.  A  non-propagating  crack  with  no 
prior  loading  history  or  existence  of  residual  strains  behind  the  crack 
tip  cannot  provide  the  closure  levels  found  in  a  propagating  crack  in  a 
finite  element  analysis.  A  plastic  wake  over  a  finite  crack  length  is 
necessary  to  develop  and  demonstrate  crack  closure. 

Under  most  loading  conditions,  the  formation  of  this  plastic  wake 
restricts  crack  surface  opening  displacements  and  reduces  the  total 
strain  occurring  immediately  in  front  of  the  crack  tip.  These  effects 
were  observed  in  the  finite  element  analysis  over  the  region  of  crack 
growth  from  the  original  crack  tip  to  its  final  location.  If  a  specimen 
has  been  'precracked'  before  the  onset  of  fatigue  loading  without 
allowing  for  the  development  of  a  plastic  wake,  then  displacements  along 
the  precrack  will  be  less  and  total  strains  observed  will  be  greater 
than  one  exhibiting  a  plastic  wake  effect.  A  'precracked'  specimen 
exhibiting  no  crack  growth  history,  either  experimentally  or 
analytically ,  will  not  accurately  demonstrate  the  effects  of  plasticity 
induced  closure  on  fatigue  crack  growth. 

The  investigation  conducted  as  part  of  this  research  effort 
indicates  that  this  plastic  wake  effect  develops  rapidly  for  fully 
reversed  loading  but  requires  longer  propagation  distances  to  develop 
under  a  stress  ratio  of  R  -  0.1.  As  a  result,  plasticity  induced 
closure  is  also  observed  to  develop  more  rapidly  at  higher  stress  levels 
than  lower  stress  levels.  Further,  the  steady  state  closure  loads,  as  a 
fraction  of  maximum  load,  decrease  with  increasing  stress  or  decreasing 
stress  ratios. 

Because  of  this  build-up  of  plastic  material,  residual  compressive 
forces  are  exhibited  during  the  unloading  portion  of  the  fatigue  cycle. 
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This  causes  an  increase  in  the  amount  of  tensile  load  required  to  open  a 
crack  and,  in  a  like  manner,  a  decrease  in  the  amount  of  external 
compression  required  to  keep  the  crack  closed.  This  plastic  wake  effect 
increases  the  values  of  P(clo)  and  P(ope),  the  closure  and  opening 
loads.  This  effect,  however,  is  not  the  same  for  all  cases.  As 
discussed  previously,  the  variations  in  crack  opening  stress  intensities 
K(ope)  with  maximum  stress  intensity  K(max)  is  neither  a  single-valued 
function  nor  a  monotonically  increasing  function  for  short  cracks  under 
different  stresses  or  stress  ratios.  For  either  value  of  R 
investigated,  K(ope)  decreases  with  increasing  stress.  Opening  and 
closing  stress  intensities,  K(ope)  and  K(clo) ,  tend  to  diverge  with 
increases  in  K(max)  or  crack  length  under  constant  amplitude  fatigue  at 
R  -  -1.0. 

The  determination  of  closure  and  opening  loads  requires 
displacement  measurements  at  some  point  along  the  crack  length.  This 
investigation  showed  that  closure  is  more  easily  determined  using 
near- field  displacements  (those  immediately  behind  the  crack  tip)  than 
far-field  (those  at  or  near  the  crack  mouth).  Differential  displacement 
plots,  however,  ease  the  determination  of  opening  and  closing  loads. 

For  high  stress  levels  where  significant  hysteresis  is  observed,  the 
inflection  point  in  the  unloading  curve  correlates  with  the  closure 
load.  These  closure  and  opening  loads,  using  either  displacement 
measurements  or  differential  displacement  plots,  do  not  depend  on 
displacement  measurement  location,  but  the  sensitivity  and  resolution 
required  to  determine  these  loads  is  greatly  affected  by  location. 
Near-field  measurements  provide  greater  sensitivity  to  closure  load 
determination. 
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These  above  observations  regarding  the  effects  of  plasticity 
induced  closure  on  fatigue  crack  growth  are  consistent  over  a  wide  range 
of  load  conditions.  High  compressive  forces,  however,  which  occur  in  a 
reversed  loading  condition  tend  to  remove  some  of  these  effects.  The 
recovery  of  some  plasticity  as  a  result  of  these  forces  along  the  crack 
surface  reduces  the  size  of  the  plastic  wake  and  reduces  the  restriction 
of  strain  related  conditions.  These  results  appear  highly  dependent  on 
the  loading  history  of  a  specimen  and  are  not  easily  predicted. 


V .  Summary  and  Conclusions 


Visco  II,  a  plane  stress/plane  strain  finite  element  code  has  been 
developed  from  the  original  Visco  program  developed  by  Hinnerichs  [71]. 
This  modified  version  incorporates  non-linear  viscoplastic ,  as  well  as, 
elastic-plastic,  material  behavior,  changing  boundary  conditions 
resulting  from  crack  growth  and  crack  closure,  and  a  data  storage/ 
restart  capability.  The  overall  solution  technique  is  the  residual 
force  method  which  maintains  a  constant  stiffness  matrix  and 
incorporates  plasticity  effects  through  the  inclusion  of  a  plastic  load 
vector.  Non-linear  material  behavior  is  handled  by  the  Bodner- Partoin 
viscoplastic  Flow  Law,  and  an  Euler  extrapolation  scheme  is  used  to 
integrate  these  relations  with  respect  to  time.  Further,  a 
Gauss-Siedel  iterative  technique  is  used  to  solve  for  equilibrium 
conditions . 

This  finite  element  code  was  employed  to  study  the  effects  of 
plasticity  induced  closure  in  short  cracks  by  analyzing  single-edge 
cracked  specimens  (a  «  .001  in)  under  various  load  conditions. 

Materials  studied  included  TI-6246  at  room  temperature  and  Inconel  718 
at  12()0°F. 

The  following  are  summary  statements  and/or  conclusions  based  upon 
the  short  crack  fatigue  analyses  discussed  herein: 

1.  A  procedure  was  developed  which  accurately  handles  crack 

closure  in  a  cyclically  loaded  specimen.  This  algorithm  is 
initiated  within  the  finite  element  code  when  negative 
displacements  occur  along  the  crack  surface.  A  vector  of 
closure  loads,  which  is  required  to  maintain  zero  displacements 
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where  closure  has  occurred,  is  then  incorporated  into  the 
equilibrium  solution  such  that  refactorization  of  the  stiffness 
matrix  is  not  required.  This  maintains  the  overall  residual 

% 

force  method  solution  technique.  This  procedure  produced 
substantial  improvements  in  the  accuracy  of  past  finite  element 
analyses  of  cyclically  loaded  center  cracked  specimens. 

2.  A  procedure  to  accomodate  a  data  storage/restart  capability 
was  incorporated  into  the  Visco  II  finite  element  code.  This 
modification  allows  for  extensive  finite  element  analyses  which 
were  previously  unattainable  due  to  computer  restraints. 

3.  The  crack  growth  algorithm  included  in  previous  forms  of  the 
Visco  finite  element  code  was  modified  and  successfully  used  to 
analyze  a  cyclically  loaded  specimen  for  the  first  time. 
Previous  analyses  used  this  algorithm  solely  to  study  creep 
crack  growth.  The  change  in  plastic  strain  exhibited 
immediately  in  front  of  the  crack  tip  on  subsequent  load  cycles 
was  used  as  a  criterion  for  crack  growth.  This  significantly 
reduced  computational  requirements  and  satisfactorily  modelled 
material  behavior  in  a  specimen  exhibiting  fatigue  crack 
growth . 

4.  Material  constants  for  TI-6246  at  room  temperature  were 
determined  for  use  in  the  Bodner  model  following  procedures 
similar  to  those  developed  by  Beaman  [72]  and  Stouffer  [ 74 ] . 
These  values  represent  a  rate- insensitive  material  at  room 
temperature  and  demonstrate  that  the  Bodner  model,  which 
incorporates  both  rate  dependent  plasticity  and  creep,  can  be 
successfully  applied  to  such  a  material.  These  constants  were 
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determined  based  upon  minimal  experimental  stress/strain  data 
and  were  subsequently  shown  to  satisfactorily  model  the 
particular  material  behavior. 

5.  The  Visco  II  code  was  used  to  model  the  growth  of  small  cracks 
to  represent  the  behavior  of  TI-6246.  The  numerical  procedure 
allowed  the  study  of  the  formation  of  a  plastic  wake,  residual 
stresses,  and  deformation  both  behind  and  ahead  of  the  crack 
tip.  This  constitutes  the  first  detailed  analytical  study  of 
the  growth  of  short  cracks  from  a  free  surface  and  the  effects 
of  plasticity  induced  closure. 

6.  At  maximum  tensile  load,  the  non-propagating  crack  is  fully 
open  and  maximum  displacements  in  the  vicinity  of  the  crack  tip 
occur  at  the  node  directly  behind  the  tip.  The  displacement 
profile  of  the  propagating  crack,  however,  shows  that  the  crack 
remains  closed  for  several  nodes  behind  the  crack  tip  and  then 
opens  gradually.  Under  cyclic  loading,  the  displacements  take 
on  unusually  high  values  directly  behind  the  crack  tip  when  no 
history  of  crack  extension  is  considered. 

7.  Plastic  strains  immediately  ahead  cf  the  crack  tip  are  higher 
for  the  non-propagating  crack  than  for  the  propagating  crack. 
Under  fully  reversed  loading,  a  severe  compression  develops  in 
the  plastic  wake  and  is  most  compressed  at  higher  stress  levels 
directly  behind  the  crack  tip. 

8.  Finite  element  modeling  of  plasticity  induced  closure  in  short 
cracks  in  fatigue  requires  some  amount  of  crack  growth  to 
develop  closure.  A  non-propagating  crack  with  no  prior  loading 
history  or  residual  strains  behind  the  crack  tip  cannot  provide 
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the  closure  levels  present  in  a  propagating  crack  using  finite 
element  simulations.  The  plastic  wake  effect  develops  rapidly 
for  fully  reversed  loading  but  requires  greater  crack  extension 
to  develop  under  a  stress  ratio  R  -  0.1.  Closure  also  develops 
more  rapidly  at  high  stress  levels  than  at  low  stress  levels. 
Steady-state  closure  loads,  as  a  fraction  of  maximum  load, 
decrease  with  increasing  stress  or  decreasing  load  ratio. 

9.  The  variation  of  crack  opening  stress  intensity,  K(ope),  with 
maximum  stress  intensity,  K(max),  is  neither  a  single-valued 
function  nor  a  monotonically  increasing  function  for  short 
cracks  under  different  stress  levels  or  stress  ratios.  For 
either  value  of  stress  ratio  considered,  R  -  -1.0  or  R  -  0.1, 
K(ope)  decreases  with  increasing  stress  as  the  nominal  applied 
stress  moves  from  0.6  to  0.9  of  the  material  yield  stress.  In 
addition,  opening  and  closing  stress  intensities,  K(ope)  and 
K(clo) .  tend  to  diverge  with  increases  in  K(max)  or  crack 
length  under  constant  amplitude  fatigue  at  R  -  -1.0.  This 
divergence  is  highest  at  a  nominal  stress  value  of  0.9  material 
yield  stress. 

10.  Closure,  as  determined  from  deviations  from  linearity  in  a 
load-displacement  plot,  is  more  easily  determined  using 
near-field  displacements  (immediately  behind  the  crack  tip) 
rather  than  from  far- field  ones  (at  the  crack  mouth)  in  small 
cracks  where  the  plastic  zone  size  is  of  the  same  order  as 
crack  length.  Differential -displacement  plots  make  it  easier 
to  determine  opening  and  closing  loads.  For  high  stress  levels 
were  significant  hysteresis  is  observed  in  the  load-displace - 
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ment  curve,  the  inflection  point  in  the  unloading  portion  of 
the  curve  correlates  with  the  closure  load  determined  from  the 
first  contact  of  a  node  behind  the  crack  tip  in  the  finite 
element  analysis.  Closure  and  opening  loads  do  not  depend  on 
displacement  measurement  location,  but  the  sensitivity  and 
resolution  necessary  to  determine  these  loads  are  greatly 
affected  by  the  measurement  location.  Near-field  measurements 
provide  this  greater  sensitivity  and  resolution  in  closure  load 
determination . 

11.  The  same  general  characteristics  regarding  the  effects  of 

plasticity  induced  closure  in  the  growth  of  short  cracks  under 
fatigue  loading  for  TI-6246  at  room  temperature  (rate 
independent  material  behavior)  were  observed  for  Inconel  718  at 
1200°F  (rate  dependent  material  behavior).  Viscoplastic 
effects  were,  however,  observed  when  considering  the  difference 
in  cyclic  load  frequency.  At  the  higher  load  frequency, 
plastic  strains  are  less  along  the  crack  length  than  for  the 
lower  cyclic  frequency.  The  size  of  the  plastic  wake  is 
smaller  at  higher  frequencies  resulting  in  less  residual 
compressive  stresses  at  the  higher  frequency.  Crack  opening 
and  crack  closing  loads  are  therefore  reduced  as  load  frequency 
is  increased. 
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Appendix  A:  Determination  of  Bodner  Constants  for  TI-6246 


The  constants  used  in  the  Bodner  model  are  both  material  and 
temperature  dependent  and  must  therefore  be  evaluated  for  each  material 
at  the  specific  temperature  of  interest.  Material  constants  for 
TI-6246  at  room  temperature  have  not  been  previously  determined  and 
must  therefore  be  evaluated  for  use  in  this  research  effort.  The 
'procedure  for  evaluating  the  Bodner  constants  developed  by  Beaman  [72], 
which  built  upon  the  works  of  Stouffer  [74],  is  used  in  conjunction  with 
experimental  stress/strain  data  provided  by  the  WRDC  Materials 
Laboratory  (Figure  2.1). 

An  initial  assumption  is  the  absence  of  thermal  recovery  of  work 
hardening.  The  governing  equation  for  determining  the  rate  of  work 
hardening, 

Z  -  m(Z1-Z)aeP  -  AZ1[(Z-Z2)/Z1]r  (A.l) 

can  be  written  as 

dZ  -  m(Z, -Z)dW  (A. 2) 

1  P 

This  equation  is  then  integrated  using  the  initial  condition  that 
Z  -  Z„  at  W  -  0  to  obtain 

0  p 


dZ/fZ, -Z)  -  mdW 
'  1  P 

(A. 3) 

ln(Z, -Z)  -  ln(Z, -Z_)  -  mW 

1  1  u  p 

(A. 4) 

Appropriate  strain  rate  data,  graphed  as  ln(Z^-Z)  versus  is  then 
linearly  approximated  to  obtain  values  of  m  (slope)  and  Zq 
(y- intercept) . 


158 


The  data  used  in  this  procedure  was  taken  from  the  experimental 
stress-strain  curve  provided.  The  elastic  modulus,  E,  is  a  reliable 
description  of  the  material's  elastic  behavior  and  values  of  stress, 
strain,  and  strain  plus  'offset'  strain  were  determined  from  the 
experimental  data  and  deviations  from  a  plot  of  the  elastic  modulus.  In 
order  to  calculate  the  plastic  work,  W^,  plastic  strains  were  first 
derived  from  the  total  strains.  This  was  accomplished  using  Hooke's  Law 
to  derive  the  plastic  strains  by 

€P  -  €  -  cr/E  (A.  5) 

Using  the  values  obtained  here,  the  plastic  work,  W^,  was  computed  using 
the  general  form 

w  n  -  w  H'1  +o-n‘1(eP(n)  -  + 

P  P 

((eP(n)  •  -r”'1))^  (a. 6) 

The  x-coordinates  of  Equation  (A. 4)  were  thereby  determined. 

Evaluation  of  the  y-coordinates  begins  with  the  evaluation  of 

Z  -a[(2n/(n^l))ln(2D0/v/3lP)]n/2  (A. 7) 

which  requires  knowledge  of  the  plastic  strain  rates.  This  was 
accomplished  using 

€e  -ff/E  -  ACT/AtE  (A.  8) 

where  €-  +  eP  and  €  «  1.3  x  10  ^  in/sec.  Since  'the  experimental 

total  strain  rate  is  held  constant,  the  time  between  data  points  can  be 
computed  using  the  change  in  total  strain.  Values  for  the  plastic 
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strain  rates  are  then  computed  over  subsequent  data  point  intervals. 

4 

In  this  procedure,  Dq  was  assumed  to  be  10  sec  and  the  values  of 
n  and  2^  were  estimated  using 

ln(-ln(v/3€P/2D0>)  -  -2nln<J  +  (2nlnZ1  +  ln((n+l)/2n)  (A. 9) 

The  linear  relationship  between  ln( -ln(v/3lP/2DQ) )  and  lnCT  permits  the 
determination  of  n,  and  subsequently  Z^.  Since  only  one  set  of 
stress-strain  data  was  available,  the  values  obtained  were  only 
approximations  of  the  actual  material  -values.  Using  this  initial  value 
of  n,  subsequent  values  were  taken  in  the  vicinity  of  this  first 
approximation  and  corresponding  values  of  Z^  evaluated.  Employing  these 
values,  subsequent  evaluations  of  Equation  A. 7  led  to  several  associated 
values  of  Z  and  ln(Z^-Z) . 

At  this  point  a  two-element  finite  element  model  employing  the 

above  estimated  values  was  used  to  evaluate  the  Bodner  constants.  Since 

it  was  assumed  that  there  was  no  thermal  recovery  present,  the  values  of 

Z^,  A,  and  r  were  set  equal  to  zero.  Additionally,  Dq  was  assumed  to  be 
U 

10  sec.  The  value  of  n  was  varied  within  a  region  near  the 
approximated  value  determined  above  and  a  resulting  value  of 
obtained.  Using  these  ' wo  values,  the  constants  m  and  Zq  were 
evaluated.  This  procedure  was  repeated  until  the  sec  of  Bodner 
constants  resulted  in  stress  -  strain  data  which  resembled  the 
experimental  data  as  nearly  as  possible.  These  resultant  constants  are 
shown  in  Table  2.2  and  the  subsequent  stress-strain  curve  in  Figure  2.1. 

The  set  of  Bodner  constants  evaluated  in  this  manner  is  only  an 
approximation  of  the  actual  material  values.  The  accurate  determination 
of  these  coefficients  requires,  as  a  minimum,  several  experimental 
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stress-strain  curves  at  high  strain  rates  with  no  thermal  recovery,  and 
several  data  points  involving  lower  strain  rates  [72],  The  values  of 
the  coefficients  determined  above  are  highly  subjective.  As  indicated, 
however,  they  produce  a  satisfactory  representation  of  the  experimenta¬ 
tion  provided. 
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Appendix  B:  Visco  II  Propram  Modifications 


The  finite  element  code,  Visco  II,  is  a  modification  of  the 
program  developed  by  Hinnerichs  [71]  in  order  to  satisfy  the 
requirements  of  this  research.  These  changes,  although  important  to 
this  study,  do  not  interfere  with  the  computational  procedures  of  the 
Bodner-Partom  Viscoplastic  Flow  Law  and  the  various  solution  techniques 
such  as  the  Residual  Force  Method  and  the  Gauss-Siedel  iterative 
algorithm.  The  modified  program  has  been  verified  by  re-analyzing 
several  cases  previously  studied  by  Henkel  [68]  and  by  comparisons  with 
works  by  Mercer  [73],  The  modifications  that  are  presented  here  are 
the  data  storage/restart  capability  and  the  crack  closure  algorithm. 

The  computer  code  involved  in  the  data  storage/restart  capability 
is  presented  in  Figure  B.l.  This  procedure  requires  the  input  of  the 
variable  NCONT  (00037)  which  indicates  if  data  will  be  stored  at  the 
completion  of  a  given  analysis  and  whether  or  not  this  particular 
analysis  is  a  continuation  of  a  previous  computer  run  (0  -  not  used, 

1  -  data  storage  only,  2  -  both  data  storage  and  program  restart 
required).  If  the  computer  analysis  is  exercising  the  restart 
capability,  the  input  data  is  read  from  storage  (00042-00083).  New 
values  of  TMAX  and  STOPCY  (00039-00041)  must  additionally  be  input  at 
this  point  to  place  time  limits  on  the  current  analysis.  If,  at  the 
end  of  the  current  run,  data  will  be  stored  for  future  use,  then  this 
is  accomplished  after  the  current  analysis  is  complete  (00371-00412). 

The  procedure  shown  here  has  been  verified  by  duplicating  Henkel's 
results  [68]  for  a  center  cracked  panel  under  load  conditions  of  - 
45  KSI-IN1/2  and  R  -  -1 .0 
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ODD  34C 

READ  AND 

00Q35C 

00036 

WRITE  Cb* 

003  37 

READC5,* 

00338 

I FC  C  NCON 

00039 

Read  c  5** 

00040 

PRINT  C£ , 

00041 

- 

PRINT (  6, 

09042 

READCl,* 

00043 

READCl,* 

00044 

READCl,* 

03045 

- 

READCl,* 

00046 

READCl,* 

00047 

READCl  ,* 

00043 

READCl,* 

00049 

READ  Cl,* 

03050 

READCl,* 

00351 

READCl,* 

00352 

READ  Cl  ,* 

03353 

READCl,* 

03054 

READ  Cl,* 

00055 

READCl ,* 

30356 

READ  Cl,* 

03057 

READCl ,* 

00059 

READCl,* 

0  00  59 

READCl,* 

0  03  63 

READCl,* 

03061 

READCl,* 

00362 

READCl,* 

0  03  63 

1  X  (.  DD  CI>, 

90064 

31  32.'  J 

00365 

READ  Cl,* 

00366 

READCl,* 

2  00  67 

823 

CONTINUE 

•’•30  63 

DD  83:  J 

0  33  69 

READCl,* 

00070 

READ  Cl,* 

0  03  71 

9  33 

CONTINUE 

030  72 

TrC I CR  »E 

03073 

Du  840  J 

30074 

READCl,* 

030  75 

34-1 

CONT INUE 

033  7s 

9  43 

READCl  ,* 

0  30  77 

IF  CSCLO 

09079 

DO  346  I 

0  30  79 

84o 

READCl,* 

000  83 

947 

CDNT INUE 

00081 

03  850  J 

00092 

READCl ,* 

003  83 

9  53 

CDNT INUE 

00084 

CALL  OUTI 

00385 

GO  TO  301 

03086 

813 

CONT INUE 

URITEC6.99) 

READC5,* I NCONT 

IPCCNC0NT.EQ.3 ).0R .(NCONT.EQ.l ) )  GC  TO  810 
R£ADC5,*)TMAX,STCPCY 
PRINTC£,*I«.NE.C  MAX  TIME  s  ",TMAX 
PRINTCfc.O^NEW  ST  CP  CYCLE  =  ••STOPCY 
READCl,*) NPATA ,M£SH£ *MESHN ,MESH3L C *SC AL*MPRI NT 
R£ADC1,*)PS,MAT,DZ2,EN,Z1,Z0,ZI  ,E**3.RN,AC 

r:adci,*inumbc  ,ncpin  »nop  in  »ncycm,  toler  ,xfac 

R£ADCl,*)VP,riELD,HP,VC*ALPHA,RAM,BERG,DTMAX,TPRINT,TPLOT 
READC1,*>NC*PMAX,PC,P£RI0D,TSTR£SS,TEPS,VH,P£,SL0PEA 
R£ADC1 , *  I NF,FR  ATE, KCC£ » OPE ,CEPT , RT ,CONV» ICRR , DTI NI T 
READ!!,*  > ICR  * J PATH *T*TP*NUMEL*NU MNP. C JINTsDDT.DT »P I »PP  *NN , TPL 
READ  (1  *  *  )  (  LM  (  I  >,1  =  1,3) 

READC1,*ICAREACII ,I=l,NUMELl 

READCl, OCNPIC  I>,NPJCI>*NPKCI)*I=1*MUMEL> 

READ  (1  ,*  ICSIGXXCII  .SISXYCI  >  *  SI  GY  Y  (  I  >  ,  SIGZZ  (  I)  ,1=1,NUMEL> 

REA  DC ! ,*>CDSIGXXCII,DSIGX Y Cl I , 031  GY YC 1 1 » DSI GZZC I  I » 1  =  1 , NUME L I 
READ  Cl,*  HDXCI  ),DY(I),DXYCI)«DZCI)«I=1*NUMEL) 

R:ADC1,*>CW»P(  II,.T>CC!>«WP~C1>,ET<I>,VUCI>,THCX>  ,I=1,KUMEL) 
READCl,*  )  C  rVPX  (I  I*  EVPr  Cl  I,  EVPXYC  1 1  ,£VPZ  Cl  1,1=1,  NUMED 
RE AD  Cl ,*  >  <E»EFF(I> ,DP£FFt 1 1 ,1  =  1 ,NUMEL> 

R  E AD  C 1 »*  I  C  NPBC  1 1  * NFI XC I ) *  SLOPE C I  I , I  =  1 ,NUMBC) 

RE  AD  Cl,* I CDSVC 1 1 ,03Y C 1 1  * X L CADC 1 1  * YLO ADC  1 1 • MAP C 1 1  *1=1 ,NUMNPI 
READC1,* I C  XORD  C 1 1  * Y  JRDC 1 1  *  DXDT Cl  I ,3 YDT Cl > *  1  =  1 *NUMNP ) 
READC1,*ICFXCI  >,FYCI),NFAC I) ,NFN ( 1 1  *  1  =  1* NFI 
READC1,*ICS<XC(I)*SXYCCII*SYXCCI I *S Y YC C 1 1 • TC R ACK ( 1 1  * 

*t.DD  CI>, 1  =  1,2?) 

dt  32:  j=:,9 

READ  Cl,*  >  CSX  VC I*J)«SXYCI,J)«SYXCI,J)*SYYCI*J)*I  =  1, MUMNP > 

REA3C1** ) (NPCI *J>*1=1* NUMNP I 

CONTINUE 

DD  837  J=1,1C 

READCl,*  > JELCJ  I 

READ  Cl,* )C JVC J ,1) ,XNCJ,I ) ,YNCJ,I 1*1  =  1*131 

CONTINUE 

IrCl CR.EQ.C)  GC  TO  845 


TO  847 


Figure  B.l 


Vises  II  -  Data  Storage/Restart 
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0  0370 
00371 
0  0372 
03373 
03374 
00375 
00376 
00377 
00378 
00379 
00380 
00381 
00382 
00383 
00384 
">385 
i  C386 
00387 
00386 
00389 
00393 
00391 
00392 
03393 
00394 
03395 
0039S 
00397 
::393 
00399 
0  3403 
00401 
00402 
03403 
00404 
00405 
0  040  i 
00407 
00403 
00409 
004  10 
00411 
03412 
00413 
00414 


IF(NCONT.£Q.C»  60  TO  890 

-  MRITEdt • ) NOATA,8£SH£tMESHN*8cSHOLO*SCALt8PRINT 
HRITE(2,*IPS,NATtD?2*EN*ZltZS*ZI*CH3»RN,AC 
MRIT£<2»*)NUM3CfNCPIN»N0PlNtNCYCMtT  OLERtXFAC 
HRITECfO  VP  *YI£LOtHPtVC»  ALPHA  *RAM*BERG»DTHAX*TPRINT*TPLOT 
WRIT£C?**>NCfPFAX*PO*P£RlOD*TSTRESS»TEPStVH#PE*SLOPEA 
HRITE(2»*JNC  *F  RATE • NODE  *OPE*CEPT #RT #C ONV  * I CRR  *DT IN IT 
«91TE<’**>i:R»JPATH*T*TPtNUH*:L*NUMNP#CJlNT*00Tt0T,PI,PP,HN,TPL 
WRITE(?«*MLM(  I>»I=1*3» 

WRIT £(?«*)< AREA! I >#I=1*HUMELI 

HRITECZt*) ( NPI (1) «NPJ(I)tNPK(I >  *I  =  1|NUHEL> 

WRIT£(2t*MSIGXXCI ) »SIGXYC I)  •SIGYV(I)«SIGZZ(I)«I=1*HUHEL) 
HRIT£<2**MDSI  GXX (  I »  *DS IGX Yl  1 1  *DSIGYY  C I)  *OSI GZZ<  I  )  *1  =1  •  NUMEL  I 
HRlTE(i»*MOK(  !ltDY(  I )  , DX Y  <  I  ) t  02  C  I )  ,I=1,NUHEL) 
URITE(2»*MUP?CII  *  HPC<  I  >  *  UP£  (I  )  >  ET C  I  MXUCI  J  «  THC I  1 «  1  =  1 » NUHE L ) 
URITEC’t-MEVOXCn  tEVP Y ( I ) tE V®XY < I > *E VPZ Cl  I  *  I =1*  NUME  L> 

WRITE C?**) (EPEFFC I ># DP EFF C I > ,1 =1 , NUHEL > 
WRITE<2i*MNP3<Il*KFIX(II*SLOPE(I)tI  =  l*NUHBC) 

HRlTECCiOCDSX  Cl  1 1 DSYC  I  )  *  XLOAD  C I  >  « YLOAD C I  >  *NAPC  I  >  »  1=1  * NUPNP ) 
WRIT£(?,*MXORO(II  ,YORDCI>  *OXDT C  I  »  » DY DTf  I  »  » I  =1  *NUHNP  > 

WRIT  £<?#•> (FX( I) *FY(I>  »NFACI)*HFNCII#I=1»NF> 

URIT£(2**MSXXC(I>  «SXYCCI>  tSYXCCI  >tSYYC<  I»  tTCRACKCI)  * 

! YL00CI)tI=l*29> 

DO  870  J=l,9 

URITE(0**MSXX  CItJI.SXYCIf  J>  tSYX  CI*J)*SYYCI*  J>*I=1  *HUMNP) 
URlTE(2t*MNP(  I«J)*I=1* NUMNP  > 

870  CONTINUE 

DO  375  J=l*10 
V.  IT"(2**)JEL(J> 

WRIT£(?»*MJN(J*I)  »XN(J»I>  «  YNC  J  »  I  )*I  =  1*13) 

375  CONTINUE 

IFCICR.EG.CI  6 C  TC  882 
DO  33.  J=l»4 

HR ITE(?»« ) (NCR  <I»J)»I  =  1 tICR> 

890  CONTINUE 

93  2  WRITE("t*)CTOL  *NCLC 

I F  ( NCIO.EQ.CI GO  TC  884 
DO  8  33  1=1  *NCLG 

983  WRITEC'tO (NCLCSC I*J)«J=lf5> 

884  CONTINUE 

DO  985  J=1»S 

HR  IT  £(*»»•  MB  (I  f  J»  »  DC  I*J)»S(  I»J)  »I=1*6> 

885  CONTINUE 
990  CONTINUE 

ENO 


Figure  B.l  (cont) 


Visco  II  -  Data  Storage/Restart 


The  crack  closure  algorithm  developed  in  this  computer  code  is 

divided  into  two  parts  and  these  are  presented  in  Figure  B.2.  The 

first  of  these  is  the  input  data  (00171-00179)  which  includes  the 

closure  tolerance  (CTOL) ,  the  number  of  nodes  that  will  be  checked  for 

closure  (NCLO) ,  and  a  matrix  (NCLOS(N.J))  identifying  the  specific 

nodes  of  interest  and  the  surrounding  elements. 

The  second  portion  of  the  closure  modification  is  the  additional 

set  of  computer  statements  within  the  equilibrium  solution.  In  lines 

00808  to  00814  the  nodes  identified  are  checked  for  negative 

displacements.  If  negative  displacements  have  occurred,  they  are  reset 

to  zero  and  resultant  closure  loads  are  calculated  (00846-00881)  . 

Convergence  of  the  closure  loads  is  checked  in  lines  00882-00885.  If 

the  difference  in  closure  loads  is  within  the  closure  tolerance,  CTOL, 

then  the  analysis  continues.  If  not,  the  program  reinitiates  the 

equilibrium  solution  and  subsequent  recalculation  of  closure  loads  based 

upon  a  new  set  of  nodal  displacements. 

As  previously  stated,  these  algorithms  have  been  verified  through 

comparisons  with  other  works.  The  results  obtained  using  the  above 

1/2 

procedures  agreed  well  with  Mercer  [73J  when  analyzing  a  35  KSI-IN 

loaded  center  cracked  panel.  Additionally,  an  analysis  of  a  45 
1/2 

KSI-IN  loaded  center  cracked  panel  agreed  with  the  results  obtained 
by  Henkel  [68],  In  both  cases,  the  use  of  the  above  algorithm 
significantly  reduced  the  residual  forces  in  the  equilibrium 
solutions . 


165 


00143 
00149 
001  SO 
00151 
09152 
00153 
00154 
00155 
99155 
00157 
00158 
00159 
00169 
00161 
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00165 
0016S 
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read  (5»oicr«tcrack(i> 

PRIST  f  6  THE  NUMBER  OF  NODES  TO  BE  CR ACKEDS"« ICR 

IF  (ICR  .EQ.  9  >G0  TO  94C 

DO  20 C C  I  =  1»ICP 

READ  (5f • > (NCR  (I«K >»K  =  1 *4) 

2000  PRINT  <6»*l(NCR(ItK)tK=l»4) 

READ  (5t*  )(  TCRACKd  >.  1  =  1  *ICR> 

DO  930  I =1 t I CR 

PRINT(6**I"N0DE  ■ *NCR (I ♦ 1 > *"POP S  AT  TINE  •tTCRACKCII 
930  CONTINUE 
943  CONTINUE 

READ  CEt*»J»ATH 

PRINT(6#*)"NUMEER  OF  J  INTEGRAL  PATHSs"* JPaTH 
IF  (JPATH  .EQ.  0I5C  TO  910 
DT  93C  1=1, JPATH 
READ  <rt*)JEL<I> 

°RINT  (6»*>"N'JR3ER  OF  ELEMENTS  IN  PATH" t I  t "EQUALS"  *JEL(I  > 
N=JEL(I» 

READ  (E«122  >CJMI«K>  *¥  N<  1 t  K)  *  Y  N  (  I  •  K  )  *  K=1  «  M  ) 

123  FORMAT  (I5.2E13.2I 
DO  925  K=1*M 

92  5  PRINT  <6«0"JN="*JMI*KI*"  XN=  "*XNCItK>*"  TN=  "*VAdtK> 

903  CONTINUE 
9 1  **  READ(:  ,•  ICTOLeKCLC 

PRINT (f»*l"CLOSURE  TCLERAN2E=  ".CTOL 
PRI NT(6t*»"N0DES  ALONG  CRACK=  "*NCLO 
IF(NCL?.Ee.?>  GO  TO  925 
03  92C  N=1 tNCLC 

r:ad<"**mnclos(n*  ji  »j=i*5» 

923  PRINT (tt*)"NGDE=  ■ tNCL9S( N*1 >*•  ELEMENTS  ARE:  "«NCL0S(N*3 >• 
1"  "*NCL0S(N,4> »"  "*NCL0S(Nf5> 

925  CONTINUE 


Figure  B.2 

Visco  II  -  Closure  Algorithm 
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DIMENSION  &(£>*£>  *D(6*6> 

1CL3S=0 

00  1  3 if  N=1  «NUHNP  - 

PLOA 0 ( i t  N >  =  D • 

PL9AD  C  2*  N ) =0 • 

PL JAD(3*N)=0« 

OOT=OT 
NPRI NT  =  2 

nurp  t=n:pi N 

N'JM3?T  =NOP  IN 

NC  VC  LEO 

RTF  =  PI  +  T  *FR4TE 

IF  {  R  TF  .LT*  1. >FTOLER=RTF»TOLER 
I f  <  R  TF  •  GEI •  i. >FTOLER=TOLER 
SURO. 

imttc=: 

03  2  9C  R=1  *NUM N° 

NU1=NAD(M» 

IF  (NUH  •  EQ«  1IPRINT  <6»O*N03E  POINT  •»«,•  IS  UNCONNECTED  ■ 
IF(SxxCR»l>*$YY(M*l>  >275*293*275 

FRX  =  XLOAD<H> . 

FRr=VL0ADCH)-PL0ADC2*N) 

03  293  L=2  »NUM 
N  =  T»n,ll 

FRX=FRx-SXXCH*  L > *QSX < N )-SK YC H* L > *DS Y <N> 
FRY=FRY-SYK<M,Ll«05K<Ul-SYY(M*L»*DSttN> 

D<rSX*<H.l  l*FRX*SKY(M,l  )«FRY**05X  (  M  > 
DY=SYX(N*l»*FRX*SYY<N»l»*FRV-OSriHI 

I=<PLOtD<l,H).ST.0.|0Y=0, 

03X(M)rDSX<H)*  XFAC*DX 
93YC  M)=DSY<9|*>XPAC*0Y 
lr  CPL^ADC  l.M)  .37.0.  >  DSYC*>=3. 

03  65  1NTT=192  *215 
IF(3S>(INTTI,5E.:,)G0  70  65 
DS*( INTT  »sC. 

I  NTT  C  =  1 

PL3AD(I, INTT>=1. 

CONT INUE 

Ic  <S M,l> >20*43*20 
IF  <SYYCM*1> >285. 33*285 
S’JM=SUN*ASS<DX/SXX(Myl  >) 

63  TO  29C 

SIH  =  SUM»ABS(DY /STY  0**1  >>  -  

63  TO  290 

SUM=SU«*ABSC0X/SXX<N*1> > «*A BS <0 Y/ SY Y  CH ,  1>  >  —  -  - 

C  3  NT INUE 


Figure  B.2  (cont) 
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008463 
O0847C 
00849 
00849 
008  50 
83851 
039  52 
00853 
80S  54 
80855  140 
00858 
8085T 
03859 
09859 
09883 
00861 
03962 
03863 
03864 
0  3865 
03965 
09967 
0  38  69 
3986=* 
00373 
09971 
0  38  72 
3:8  73 
03974 
03975 
099*^5  159 
00977 
33373  155 
0  33  79 
03889  179 
00301  175 
00882  153 
00883 
308  94 
00895 
03985  195 
09997  196 
03893 
3  "9  99 
00993  185 
00991 
03892  293 
09893 
02894 
00895  205 
00896  210 
03897 
00898 


DETERMINATION  OF  CLOSED  NODAL  LOADS 

I  TER  =0 
I CLOS-1 

DO  195  M=1»NUMNP 

IF  iPLOADC  1**0  .EQ.O.  160  TO  195 

PLOAD<3*H>=0. 

DO  140  N=1«NCL0 

IP  < NCL0SCN*1 ) .EQ.N*L=N 

CONTINUE 

LN=NCL0SCL*2)*2 

00  150  N=3*LN 

I=NPI  <\:iOS(Lt  N>  > 

J=N°J(NCLOS(L»ft)> 

K  =  N®K(NCLOS(L*A)> 

CALL  BHATRIX  CMCLCS  CLtN  I  *B  1 

E’X=edtl>*DSX  tn*Bdt3)0$X<J>*BU.51*DSX<K> 
EPY=8C2*  2)*DST  (I)«9(2»4)*DST(J)4-BC2»6)*DSY(K) 
GAM=B<3.1»«0SX  ai*B(3#2l*DSYCI  > *B f 3 t 3 > *DSX C J ) 

1-*B<  3  *4  )*DSY( J) ♦  BI  3*5l*DSX(K)*Bt 3«6)*DSY<K> 

CALL  DMATRIXCNCLOSIL*NI »D> 

EPXX=EPX-EV»X< NCLCS(L*NI> 

EPYY=EPY-EVPY(  NCLC3( LtN ) ) 

SAMMsGAM-EVPXY (NCLOS(L»N> ) 

X=D(lil >*EPXX»CC1*2)*EPYY 
Y=0C2*1I*EPXX*DC2*2I»EPYY 
XY  =  D  (3»  3)*  GANN 
IFd.EQ.MIGO  TC  163 
IFCJ.EQ.MIGO  TC  165 
IF(K.E0.M1G0  TC  170 

PL0AD(3*M»=PL0ADl3#P>»3Cl*21*X»Bt2*2»*Y*BC3.2>*XY 
GO  TO  175 

?L3AD<3#M)=PL0ADC3fP)*BClt4)*X*B{2»4)*Y*B(3*4>»XYi 
GO  TO  17  5 

PL3AD(3»M)=*L3AD(3*P>*B(lt61*X*Bf2*  6>* Y^B ( 3*  6)*X Y 

CONT INUE 

CONTINUE 

JF<t»L0AD<2«M).  EQ.9.»PLOAD(2tN)=.01 

DIFF=ABS( <PLOAOC3*H»-PLOAD(2*N) > /PL0ADC2*M>> 

IF(0IFF,GT«CT0L»ITER=1 

CONT INUE 

CONTINUE 

IFdTER.EQ.3  >GC  TC  200 
DO  185  M=ltNUNNP 
PLOAD(2iMl=PLDADC3#M> 

63  TO  105 

DO  210  H=1»NUMNP 

IF  <  PLOADC  1»N>  ,E0.9.  )G0  TO  295 

PR I  NT  (  6*  O  H*  •  ■*PL0AD<3iN) 

CONT INUE 
CONTINUE 

RETURN  - 

END 

Figure  B.2  (cont) 
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